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ABSTRACT 

We construct a catalog of radio sources detected by GB6 (6 cm), FIRST and NVSS (20 cm), 
and WENSS (92 cm) radio surveys, and the SDSS optical survey. The 2.7 million entries in the 
publicly-available master catalog are comprised of the closest three FIRST to NVSS matches 
(within 30") and vice- versa, and unmatched sources from each survey. Entries are supplemented 
by data from the other radio and optical surveys, where available. All objects with even a small 
probability of physical association are included, such that catalog users can easily implement 
their own selection criteria for data analysis. We perform data analysis in the ~ 3000 deg 2 region 
of sky where the surveys overlap, which contains 140,000 NVSS-FIRST sources, of which 64,000 
are detected by WENSS and 12,000 by GB6. About one third of each sample is detected by 
SDSS. An automated classification method based on 20 cm fluxes defines three radio morphol- 
ogy classes: complex, resolved, and compact. Radio color-magnitude-morphology diagrams for 
these classes show structure suggestive of strong underlying physical correlations. Complex and 
resolved sources tend to have a steep spectral slope (a ~ —0.8) that is nearly constant from 6 
to 92 cm, while the compact class (unresolved on ~ 5" scale by FIRST) contains a significant 
number of flat-spectrum (a ~ 0) sources. In the optically-detected sample, quasars dominate 
the flat-spectrum compact sources while steep-spectrum and resolved objects contain substan- 
tial numbers of both quasars and galaxies. Differential radio counts of quasars and galaxies are 
similar at bright flux levels (> 100 mJy at 20 cm), while at fainter levels the quasar counts are 
significantly reduced below galaxy counts. The optically-undetected sample is strongly biased 
toward steep-spectrum sources. In samples of quasars and galaxies with SDSS spectra (2,885 and 
1,288 respectively), we find that radio properties such as spectral slope, morphology, and radio 
loudness are correlated with optical color and luminosity. 

Subject headings: catalogs — galaxies: active — radio continuum: general 



1. INTRODUCTION 

Quasars and powerful radio galaxies dominate 
the observed counts of continuum radio sources 
above milliJansky flux levels, and display spec- 
tacular morphological variety that is correlated 
with other properties such as spectral slope and 
luminosity. The unification paradigm for ra- 



dio galaxies and quas ars ( Urrv fc Padovani 19951: 



Jackson fc Wa ll 1999) attempts to explain much 
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of this rich variety of observational data as arising 
from essentially the same anisotropic processes 
which appear very different to us because of vary- 
ing viewing angles to the radio jets. This con- 
jecture has fundamental implications for our un- 
derstanding of quasars and galaxies, but for a 
strong test of the unification paradigm one essen- 
tially needs a large statistical sample with well- 
controlled selection criteria and robust estimates 
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of the source morphology, as well as appropriate 
models to interpret the data. 

Statistical studies of radio emission from ex- 
tragalactic sources are entering a new era, result- 
ing from the availability of large sky radio surveys 
that are sensitive to milliJansky flux levels (e.g., 
Becker et al1ll995HCondon et al.l[l99"i iDe BreucM 



2000). The catalogs based on these surveys con- 
tain large numbers of sources, have high complete- 
ness and low contamination, and are available in 
digital form. The wide wavelength region spanned 
by these surveys, from 6 cm for GB6 to 92 cm for 
WENSS, and detailed morphological information 
at 20 cm provided by FIRST and NVSS, allow sig- 
nificant quantitative and qualitative advances in 
studies of radio sources. In addition, the optical 
catalog obt ained by the Sl oan Digital Sky Sur- 
vey (SDSS: lYork et al1l2000h can be used to sepa- 
rate quasars from galaxies, and the redshifts mea- 
sured by SDSS allow a comprehensive study of the 
optical-radio correlation for quasars and galaxies. 

FIRST and NVSS, conducted separately at 
the Very Large Array (VLA), were the first ra- 
dio surveys with sufficiently high angular resolu- 
tion to allow unambiguous matching with deep 
optical surveys, providing identifications for a 
large number of radio sources. The two sur- 
veys have the same radio frequency, but FIRST 
goes slightly deeper with higher resolution and 
smalle r sky coverage ( jQl) 
(|l999h and I Sadler etall (11999 



Machalski fc Condon 
2002) measured 



the radio luminosity function of radio-loud active 
galactic nuclei (AGN) and star-forming galax- 
ies by cross-correlatin g the NVSS with spectro - 
scopic galaxy surveys. iMagliocchetti et all (2002) 
matched FIRST to the 2 degree F ield Galaxy 
Redshift Survey (jColless et alj 120011 ). using spec- 
tra to classify galaxies as "classical" radio galax- 
ies, starb urst and late - type galaxies, and Seyfert 
galaxies. Ilvezic et al. ( 2002 ) cross-correlated the 
FIRST survey with the SDSS photometric survey, 
resulting in a much larger sample of optical iden- 
tifications (about one third of FIRST sources are 
matched to an SDSS source, and 0.16% of SDSS 
sources are matched to a FIRST source, with a 
matching radius of 2"), with galaxies outnum- 
bering quasars 5:1, and a negligible frac tion of 
radio stars in the sample. iBest et al. I (|2005h cross- 
correlated the SDSS spectroscopic galaxies sample 
with both FIRST and NVSS, and found that AGN 



dominate radio counts d own to 5mJy at 20 cm . 
The more recent work of Mauch &: Sadler! lj2007h 
includes the largest radio-selected galaxy sample 
available from a single radio survey, combining 
the NVSS with th e 6 degree Field Galaxy Survey 
( Jones et al. 120041) ; ;hey confirm that radio- loud 



AGN and star-forming galaxies have quite differ- 
ent distributions in the plane of radio power ver- 
sus absolute K band (infrared) magnitude. The 
study of radio a nd optical properties of quasars 
was extended bv lde Vries et al.1 ( 20061 ) who found 
that 10% of SDSS quasars have detectable radio 
cores at 1.4 GHz (> 0.75 mJy), and 1.7% have 
double-lobed morphology, i.e., are associated with 
multiple FIRST components. 

In this paper, we describe the construction of 
a unified catalog of radio sources detected in the 
6 cm GB6, 20 cm FIRST and NVSS, and 92 cm 
WENSS radio surveys, and the SDSS optical sur- 
vey (Data Release 6, hereafter DR6). We began 
by merging the FIRST and NVSS surveys into 
a single catalog containing over 2 million sources 
detected by at least one survey; about 500,000 
sources are detected by both FIRST and NVSS. 
Where available, we supplement this survey with 
data obtained by the GB6 and WENSS surveys 
that enable the computation of radio spectral 
slopes; about 30,000 sources are detected by all 
four radio surveys. The radio sources were also 
cross-correlated with the optical SDSS catalogs; 
nearly 92,000 FIRST sources have an SDSS coun- 
terpart within 2". 

The radio and optical surveys chosen for in- 
clusion in the unified catalog primarily cover the 
northern celestial hemisphere. We have opted to 
use these surveys in order to take advantage of the 
high astrometric accuracy of FIRST, which, de- 
signed to cover the same region of sky as the SDSS, 
is limited to the northern Galactic cap. Several 
recent large-area radio surveys covering a wide 
range of frequencies are also available in the south. 
These include the Sydney University Molongl o 
Sky Survey at 36 cm (SUMSS: lMauch et al.l[2003l) . 
the Parkes-MIT-NR AO survey at 6 cm (PMN; 
Gregory et al. 1994h . and the Aus tralia Telescope 



20-G Hz (1.5 cm) survey (AT20G; |Massardi et al. 



2008) 



The main advantages of the unified catalog pre- 
sented in this paper are the multi-wavelength ra- 
dio data (92 cm and 6 cm in addition to 20 cm), 



2 



and the increased number of optical identifica- 
tions from the SDSS DR6. By combining these 
five extensive surveys, we have assembled precise 
astrometric measurements, flux at multiple wave- 
lengths, spectral index, morphological information 
(size, shape, and orientation of resolved objects), 
and optical identifications into a single compre- 
hensive catalog of radio objects. At optical wave- 
lengths, the SDSS provides colors and classifica- 
tion (i.e., into quasars and galaxies) for the ~ 30% 
of sources detected optically. The unified cata- 
log is primarily a resource of low-redshift (z < 2) 
quasars and radio galaxies with AGN, but may 
also prove to be an important source of rarer ob- 
jects such as radio stars, high-redshift quasars, and 
high-redshift galaxies. The area observed by all 
five surveys is nearly 3, 000 deg 2 . The unified cat- 
alog provides comprehensive multi- wavelength ob- 
servations at greater depth and for a larger number 
of sources than any previously available catalog of 
radio objects. 

The limits in sky coverage of the catalog are 
defined by the FIRST and NVSS sky coverage 
(Fig. []}. The unified catalog thus includes sources 
in the Galactic plane which are imaged by the 
NVSS. Galactic sources must be studied with care: 
due to the nature of interferometry, radio surveys 
(depending on their angular resolution) generally 
do a poor job of imaging highly-extended sources 
such as Galactic HII regions and supernova rem- 
nants, whose angular sizes often reach several ar- 
cmin or more. The analysis presented in this paper 
is limited to the sky covered by FIRST, which is 
greater than 30° from the Galactic plane. 

In this paper, we discuss scientific applications 
of the unified radio catalog and present a pre- 
liminary data analysis. In a companion paper 
(A. Kimball et al. in preparation), we will ex- 
pand upon and refine this analysis by comparing 
the source distribution in radio morphology, radio 
color, and optical classification sp ace to the mod- 
els of iBarai fc Wiital (|2006l . l2007h . 

The remainder of the paper is laid out as fol- 
lows. In we describe the surveys used to create 
the unified radio catalog. In ^3]we discuss the cre- 
ation of the catalog, including completeness and 
efficiency of the matching algorithms. In |2] we 
present a preliminary analysis of the radio source 
distribution according to radio morphology, radio 
color, optical identification, and, for sources with 



spectra, redshift and luminosity. We also discuss 
some catalog applications. We summarize our re- 
sults, and discuss the suitability of the catalog for 
comparison with radio evolution models, in 

2. SOURCE SURVEYS 

Before proceeding with a description of merging 
procedure for the FIRST and NVSS catalogs, we 
briefly describe each survey used in creating the 
multi-wavelength radio catalog, including the sky 
coverage, wavelength, astrometric accuracy, and 
flux limit (see Table[T]for summary). Sky coverage 
of each survey is shown in Figure [U The region of 
sky observed by all of the contributing surveys is 
indicated by the cyan solid line. 

2.1. Radio surveys 

2.1.1. FIRST 

The FIRST survey (Faint Image s of the Ra 



Becker et al 



dio Sky at Twenty centimeters, 
19951 ) used the VLA to observe the sky at 20 cm 
(1.4 GHz) with a beam size of 5"4 and an rms 
sensitivity of about 0.15mJy beam -1 . Designed 
to cover the same region of the sky as the SDSS, 
FIRST observed 9, 000 deg 2 at the north Galactic 
cap and a smaller ~ 2.5° wide strip along the Ce- 
lestial Equator. It is 95% complete to 2mJy and 
80% complete to the survey limit of 1 mJy. The 
survey contains over 800,000 unique sources, with 
astrometric uncertainty of < I". FIRST includes 
two measures of 20 cm continuum flux density: the 
peak value, -Fpeak, and the integrated flux density, 
i^int , measured by fitting a two-dimensional Gaus- 
sian to co- added images of each source. 

2.1.2. NVSS 



The NVSS (NRAO-VLA Sky Survey: ICondon et al 



19981 ) was also carried out using the VLA radio 
telescope to observe the sky at 20cm (1.4 GHz), 
the same wavelength as FIRST. However, the 
NVSS used a different antenna configuration, re- 
sulting in a lower spatial resolution (45" beam -1 ). 
Lower resolution radio surveys provide more ac- 
curate flux measurements for extended sources, 
where high-resolution surveys can miss a signifi- 
cant fraction of the flux. The astrometric accuracy 
ranges from 1" for the brightest NVSS detections 
to about 7" for the faintest detections. The survey 
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covers the entire sky north of 8 = —40° and con- 
tains over 1.8 million unique detections brighter 
than 2.5 mjy. 

For the analysis of this paper, we adopt in- 
tegrated NVSS flux densities rather than the 
peak flux densities report ed in the NVSS cata- 



log. ICondon et al.l (|1998f ) provides formulas for 



converting peak flux densities to integrated flux 
densities, as well as formulas for errors in NVSS 
measured and calculated values. The unified cata- 
log includes the integrated NVSS flux density, the 
corrected peak flux density, and the deconvolved 
major and minor axis sizes. 

2.1.3. WENSS 

WENSS (Westerb ork Northern Sky Survey; 
Rengelink et alJl997t l is a 92 cm (325 MHz) survey 
completed with the Westerbork Synthesis Radio 
Telescope. It maps the radio sky north of 5 = 29° 
to a limiting flux of 18 mjy, with a beam size of 
54" x 54"cosec(5). The complete survey contains 
almost 220,000 sources, with a positional uncer- 
tainty of < 1.5" for bright sources and < 5" for 
faint sources. 

2.1.1 GB6 

The GB6 survey at 4.85 G Hz (Green Bank 6 cm 
survey; iGreeorv et aT1ll996h was executed with the 
(now defunct) 91m Green Bank telescope in 1986 
November and 1987 October. Data from both 
epochs were assembled into a survey covering the 
0° < 8 < 75° sky down to a limiting flux of 
18mJy, with 3.5' resolution. GB6 contains over 
75,000 sources, and has a positional uncertainty 
of about 10" at the bright end and about 50" for 
faint sources. 

2.2. Optical Surveys: SDSS 

We use photometric and spectroscopic coverage 
from the sixth data release (D R6) of the Sloan 
Digital Sky Survey (SDSS: see lYork et alJ 12000b 



Stoughton et al.l l2002t lAdelman-McCarthv etal 



2007, and references therein). The survey, not 
yet finished but nearing completion, will eventu- 
ally cover 10, 000 deg 2 in the northern galactic cap 
and a smaller region on the celestial equator. DR6 
covers roughly 9, 600 deg 2 , and contains photomet- 
ric observations for 287 million unique objects, as 
well as spectra for more than 1 million sources. 



2.2.1. The SDSS Photometric Survey 

The SDSS photometric survey contains flux 
densities of detected objects measured nearly 
simultaneously in five wavelength bands (it, g, 
r, i, and z: Fukugita et al. I Il996h with effec- 
tive wa velengths of 3551 4686, 6165, 7481, and 
8931A (jGunn et al.lll998f l. The catalog is 95% 
complete for point sources to limiting AB mag- 
nitudes of 22.0, 22.2, 22.2, 21.3, and 20.5 respec- 
tively. Typical seeing is about 1.4" and positional 
uncertainty is less than about 0.1" (r ms per co- 
ordinate for sources with r < 20.5, IPier et al 



2003). The photometry is repeatable to 0.02 mag 
(rms for sources n ot limited by photon statistics, 
Ivezic et al.l [2003) and w ith a zeropoint un cer- 
tainty of ~0. 02-0.03 mag ((ivezic et alj|2004al) . A 
compendium of other technical details about SDSS 
can be found on the SDSS web siteQ which also 
provides an interface for public data access. 

Optical magnitudes in this paper refer to SDSS 
model magnitudes, measured using a weighting 
function d etermined from the o bject's r band 
image (see IStoughton et all 12002! ) . The weight- 
ing function represents the better fitting of an 
exponential profile and a de Vaucouleurs pro- 
file. The chosen model is then used to deter- 
mine the magnitude in all five bands. We cor- 
rected all magnitudes for G alactic extinction fol- 
lowing ISchlegel et al.l (|1998l ) . When selecting can- 
didate matches from the SDSS, we required unique 
sources that are brighter than r = 22.2 or brighter 
than z = 21.2. 

The morphological information from SDSS im- 
ages allows reliable star-galaxy separation to r ~ 
21.5 (jLuoton et alj|2002llscrantoii et alj|2002h . In 



brief, sources are classified as resolved or unre- 
solved according to a measure of light concentra- 
tion that represents how well the flux distribution 



matc hes that of a point source (jStoughton et al 
2002). As sources that emit strongly in the radio 



are almost exclusively extra-galactic at the high 
latitudes observed by the SDSS > 30), this 
classification effectively divides radio sources into 
"galaxies" (resolved) and "quasars" (unresolved). 
The optically-unresolved sources may also include 
a small number of galaxies that are unresolved in 
the SDSS images, as well as a small fraction of ra- 
dio stars. For the remainder of this paper, we refer 



1 http://www.sdss. org 
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to the two sets of photometric optical sources (as 
opposed to those with available spectra, see below) 
as "galaxies" and "optical point sources" . 

2.2.2. The SDSS Spectroscopic Survey 

A subset of photometric sources are chosen for 
spectroscopic observation according to the SDSS's 
spectral target selection algorithms. Targeted ex- 
tragalactic sources include th e flux-limited "main " 



galaxy sample (r < 17.77; IStrauss et a l. 2002), 
the lu minous red gal axy sample (Eisenstein et al.1 
l200lh . and quasars ([Richards et all l2002h . DR6 
contains spectra for about 100,000 quasars and 
790,000 galaxies. The spectral wavelength range is 
3800 - 9200A, with a resolution of 1800 and red- 
shift accuracy of 30kms _1 (estimated from the 
main galaxy sample) . 

The main galaxy sample includes nearly all 
(~ 99%) galaxies brighter than r < 17.77, result- 
ing in a sky density of ss 90deg -2 . Some targets 
are rejected on the basis of low surface bright- 
ness (where redshifts become unreliable or targets 
are spurious) or high flux (which can contaminate 
neighboring fibers). Due to the physical thickness 
of the spectral fibers, two galaxies closer than 55" 
cannot be observed at the same time, although 
overlapping spectral plates allow some initially- 
skipped galaxies to be picked up later. A sec- 
ond galaxy targ et algorithm selects lu minous red 
galaxies fLRGs; lEisenstein et al. I l200lh . which are 
typically the brightest members of galaxy clusters. 
LRG targets are selected primarily by color, based 
on known LRG spectra at different redshifts. The 
resulting sample is typically more luminous and 
much redder than the main galaxy sample, and 
extends to fainter apparent magnitudes than the 
main sample. 



The quasar target selection algorithm ([Richards et al 



2002) selects targets from unresolved objects with 
i < 19.1 and colors similar to redshift < 3 quasars, 
unresolved objects with i < 20.2 and colors sim- 
ilar to higher redshift quasars, and unresolved 
sources within 2" of a FIRST source (i < 19.1). 
The completeness of the sample is w 95% and 
the selection efficiency is « 66%. Contamination 
is much higher than for the galaxy sample: the 
latter consists of resolved sources, which are thus 
nearly always galaxies, whereas the quasar sample 
is contaminated by non-quasar point sources, such 
as distant galaxies or Galactic stars with quasar- 



like colors. In fact, nearly 30% of the targets turn 
out to be stars. 

To distinguish the spectroscopic SDSS sample 
from the photometric SDSS sample, we refer to 
sources from the former as "spectroscopic galax- 
ies" and "spectroscopic quasars" . 

3. THE UNIFIED CATALOG 

The cross-identification of radio surveys at dif- 
ferent wavelengths, and with d ifferent resolutions , 
is not a straightforward task dBecker et al. 1995; 
l lvezic et all 120021 : Ide Vries et all I2006L iLu et al 
l2007h . However, the accurate astrometry of 
FIRST allows simple positional matching to the 
other radio surveys, and to the optical SDSS, with 
high completeness and low contamination. Cata- 
log entries are defined by a detection in at least 
one of the two 20 cm surveys (FIRST and NVSS). 
Where available, the 20 cm data is supplemented 
with 92 cm and 6 cm radio data, from WENSS 
and GB6, and with optical observations from the 
SDSS. 

In the remainder of this section, we describe 
the matching technique used to create the unified 
radio catalog. The complicated procedure out- 
lined herein is intended to include all objects with 
even a small chance of being physically-real asso- 
ciations. Therefore, catalog users are not limited 
to the matching techniques or matching distances 
used by the authors in the subsequent analysis in 
this paper. The generous matching radii used to 
create the catalog can easily be restricted by fu- 
ture users using the catalog parameters detailed 
in Appendix [Al 

3.1. Defining catalog entries from the 
FIRST and NVSS surveys 

FIRST and NVSS both observed the sky at 
20 cm (1400 MHz). Of the four radio surveys 
incorporated into the unified catalog, these two 
have the faintest flux limits and the highest spa- 
tial resolution, with FIRST having even higher 
resolution than NVSS. Radio surveys necessar- 
ily involve a trade-off between high resolution, 
which allows for accurate determination of posi- 
tions, and low resolution, which allows for the de- 
tection of low surface brightness sources and com- 
plete flux measurements for extended sources. By 
combining high-resolution FIRST with the lower- 
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resolution NVSS, it is possible to achieve the best 
of both options. FIRST provides accurate posi- 
tional measurements, and the NVSS provides ac- 
curate flux measurements for extended and low- 
surface brightness sou rces, where FIRST underes 



timates the radio flux (Becker et al.lll995; Lu etal 



20071 ) 



3.2. Matching algorithm for the radio sur- 
veys 

Due to different spatial resolution and faint lim- 
its, combined with the complex morphology of 
radio sources, the positional matching of FIRST 
and NVSS sources cannot be done in an a pri- 
ori correct way. For this reason, we adopt a 
flexible method that enables later sample refine- 
ment after the initial positional association. The 
unified catalog contains the three closest FIRST 
matches to an NVSS source within 30", and the 
three closest NVSS matches to a FIRST source 
within 30". For each source, we also record the 
total number of matches found within 5", 10", 
30", and 120", which will allow users to inves- 
tigate radio properties as a function of environ- 
ment density. The complete catalog contains over 
2.7 million entries, including FIRST-NVSS asso- 
ciations, isolated^ FIRST or NVSS sources, and 
NVSS sources that lie outside the FIRST survey 
coverage. Because we match FIRST to NVSS and 
then NVSS to FIRST, there are necessarily du- 
plicate catalog entries: for example, a very close 
match will appear once as a FIRST to NVSS 
match and once as a NVSS to FIRST match. How- 
ever, the catalog includes parameters to easily ex- 
tract specific data samples, including the elimina- 
tion of any duplicates. Similarly, the processing 
flags can be used to treat cases such as distinct 
FIRST sources matched to the same NVSS source 
(and vice-versa, although the latter case is rare 
due to the lower spatial resolution and brighter 
limit of NVSS). For further details, we refer the 
reader to Appendix [A] 

Measured FIRST and NVSS sky positions are 
rarely exact even for the same source. As shown 
in Figure [21 FIRST has more accurate astrometric 
measurements than NVSS, due to its higher spa- 
tial resolution. We therefore use FIRST to desig- 



nate the sky position of each catalog entry, when 
possible. For NVSS sources without a FIRST 
match, we retain the NVSS coordinates. The des- 
ignated position is then used when searching for 
counterparts in the other surveys. 

To match the FIRST and NVSS surveys, we 
use a matching radius of 30". We use a larger ra- 
dius of 120" to match to WENSS and GB6, since 
they both have lower positional accuracies. We 
also record the total number of WENSS and GB6 
matches found within 120". All of the match- 
ing radii used to create the catalog are generously 
large, to ensure that the majority of physically- 
real matches are included. To select smaller, 
cleaner samples for analysis, the catalog user can 
choose matches based on distance, effectively ap- 
plying a smaller matching radius. The contamina- 
tion and completeness as a function of matching 
radius are discussed in H3.4I 

The surveys used in the creation of the unified 
catalog are themselves catalogs of individual radio 
components. Thus, large multi-component physi- 
cal sources can be resolved into separate detec- 
tions in the high angular-resolution surveys, par- 
ticularly in FIRST (5" resolution), but also in 
NVSS and WENSS (~ 50" resolution). iBest et al 



2 Isolated refers here to sources with no neighbors from the 
other 20 cm survey within 30" . 



(|2005l hereafter B05) developed a clean sample of 
double-lobed radio-loud galaxies by matching ra- 
dio components to spectroscopic SDSS galaxies. 
B05 used the optical core position and the opti- 
cal galaxy alignment to eliminate unlikely lobe- 
configurations, based on lobe opening angles and 
distances. The complexity of a matching algo- 
rithm based only on radio components, as pre- 
sented in this paper, is necessarily much greater 
because the optical position of the core is unknown 
for the majority of sources. We therefore settle 
on the matching scheme outlined above, which al- 
lows sample selection with the user's own preferred 
matching criteria, without the need to repeat the 
work that has gone into developing the unified cat- 
alog presented here. 

3.3. Matching algorithm for the SDSS 

We correlate the matched radio sources with 
the SDSS photometric survey, including spectro- 
scopic data when available, using a matching ra- 
dius of 60". In addition to the nearest neigh- 
bor, the catalog includes the brightest neighbor 
within a pre-defined radius (3", 10", 30", or 60"; 
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see Appendix [X] for details). The total number of 
matches within the same pre-defined radius is also 
recorded, indicating optical source sky density at 
each entry's location. 

We note that some lobe-dominated radio 
sources may be missed when positionally matching 
SDSS and FIRST source positions. While each ra- 
dio lobe will be included individually in the radio 
sample, it will be excluded from the optically- 
matched subset for lobe-core distances larger than 
the SDSS matching radius. Lu et al. ( 2007h esti- 
mated that about 8.1% of radio quasars do not 
show a radio core within 2" of their optical posi- 
tion, and thus would b e miss ed by this algorithm. 
Howev er, llvezic et al. I diool and Ide Vries et al.l 
( 20061 ) found fewer than 5% and 2%, respectively, 
of quasars are double-lobed in FIRST, using opti- 
cal core positions without assuming radio emis- 
sion in the core. Because NVSS has a much 
larger spatial resolution than FIRST, even fewer 
double-lobed sources can be expected in NVSS. 
B05 found 6% of their spectroscopic SDSS galax- 
ies to be double-lobed in NVSS. Note that while 
the percentage is larger for galaxies, as expected 
from orientation effects, this value is an upper 
limit for radio-optical samples using photometric 
optical data, as the spectroscopic SDSS sample 
is comprised of the nearest sources, which thus 
have large angular size. Since the effect is not 
overwhelming, the unified catalog only includes 
direct positional matches between SDSS and ra- 
dio (FIRST or NVSS) catalogs. However, catalog 
users can easily repeat the procedure s for finding 



doubl e- lobed sources develo p ed bv llvezic et al.. 
(|2002h : Ide Vries et all (|2006h : iBest et al.l (j2005h . 



which would involve matching external optical 
samples to the radio catalog sources. 

3.4. Completeness and efficiency of matched 
samples 

A side effect of using large matching radii is 
increased contamination by coincidental "line-of- 
sight" matches to physically unrelated objects. 
Using the nearest-neighbor distributions, we es- 
timate efficiency (fraction of matches which are 
physically real) and completeness (fraction of real 
matches that were found) as a function of match- 
ing radius. Figure [3] shows the distribution of 
distances between FIRST sources and the near- 
est neighbor from the NVSS, WENSS, GB6, and 



SDSS (photometric) surveys. For comparison, it 
also shows the estimated level of background con- 
tamination. The nearest-neighbor distributions 
peak at small distances, where associations are 
likely to be real, then decrease sharply; the width 
of the peaks depends on the astrometric accu- 
racy of the corresponding surveys. The distribu- 
tions slowly rise again at large distances due to 
increased background contamination. Efficiency 
and completeness are estimated using a model fit 
to the nearest neighbor distributions. Estimated 
values of completeness and efficiency based on the 
model fitting are listed in TableO For the samples 
discussed here, typical values are > 90% complete- 
ness and > 80% efficiency. 

As discussed above, the unified catalog in- 
cludes individual radio components, and does not 
correctly handle NVSS multi-component sources. 
Two lobes from a very large (> 100") source may 
appear as two detections in NVSS and thus in the 
unified catalog as well. The values in Table[2]show 
that when matching up individual radio compo- 
nents between different surveys, the completeness 
is well over 90%. As stated in §3.3[ investigations 
using optical and radio surveys suggest that only 
a few per cent of radio sources are double-lobed in 
FIRST; much smaller percentages are expected in 
the NVSS. Indeed, B05 note that the NVSS reso- 
lution is large enough that ~ 99% of radio sources 
are contained in a single NVSS component. Al- 
though the missing double-lobed sources are not 
addressed further in this paper, we tackle the issue 
in the companion paper, where we thoroughly dis- 
cuss the classification of a WENSS-NVSS-FIRST 
subsample (see §5.4.11 for discussion). It will in- 
clude the missing double-lobed sources, found by 
determining the complete NVSS-FIRST environ- 
ment around each WENSS object, requiring ex- 
tensive visual analysis. 

3.5. FIRST matching statistics 

Of the four radio surveys used to create the uni- 
fied catalog, FIRST extends to the faintest flux 
limit. We discuss here the likelihood of finding a 
FIRST source in one or more of the other three ra- 
dio surveys. In this section and for the remainder 
of this paper, we limit our analysis to the common 
region observed by all of the contributing surveys 
(see Fig.Q}. 

The number of matches found between different 
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surveys is a strong function of source flux. The 
flux range of objects in the unified catalog covers 
several orders of magnitude. Thus, for convenience 
we conve rt radio flux to an " AB radio magnitude" , 
following llvezic et al.l (|2002T ): 



iradio = -2.5 log 



( fu 



V3631Jy 



(1) 



where /j n t is the integrated flux density. This for- 
mula places the radi o magn itudes on the AB„ sys- 
tem of loke" fc Gunnl (1983). An advantage of that 
system is that the zero point (3631 Jy) does not de- 
pend on wavelength and thus enables convenient 
data comparison over a large wavelength range. 
For example, a source with a radio flux of 1 mJy 
has tradio = 16.4; if it has constant F v (a flat spec- 
trum), its visual magnitude is also V — 16.4. 

Figure |4] shows the magnitude distribution for 
FIRST sources with and without radio matches 
from the other surveys. Each panel roughly indi- 
cates the £first magnitude corresponding to the 
other survey's faint limit. For example, the frac- 
tion of FIRST sources with a counterpart in NVSS 
is greater than 0.9 at ^first = 15 but drops to 

0. 5 at £first = 15.5 as shown in Panel A, in- 
dicating that the NVSS faint limit corresponds 
to ^first ~ 15. The sensitivity limit of GB6, 
£gb6 = 13.3, occurs at £first ~ 12 (panel B). The 
limit is brighter in FIRST than GB6 because, as is 
well known and as we confirm with a much larger 
sample in N4.21 most radio sources are intrinsi- 
cally fainter at 6 cm than at 20 cm (in terms of F v ; 

1. e. corresponding faint limits depend on spectral 
slope) . Requiring a GB6 detection therefore biases 
a sample to the brightest sources. Panel C shows 
that the WENSS faint limit, £ WE nss = 13 - 3 > cor- 
responds to a limit of £first ~ 14, because radio 
sources are, again on average, intrinsically brighter 
at 92 cm than at 20 cm. 

Roughly 60% of FIRST sources have an NVSS 
counterpart within 30". Of those radio sources ob- 
served by both surveys (» 48deg -2 ), 50% have a 
WENSS counterpart and 12% have a GB6 coun- 
terpart within 120", while 11% have both. The 
sky density of sources observed by all four radio 
surveys is « 5.3 deg -2 , or ~ 5% of the FIRST 
source density. 



3.5.1. NVSS sources without a FIRST counter- 
part 

Because FIRST goes fainter than the NVSS by 
a factor of 2.5, there are many FIRST sources 
without an NVSS counterpart (~ 38deg -2 , see 
Panel A of Fig. H - 0.53 deg" 2 for FIRST sources 
brighter than lOmJy). However, the catalog also 
contains some NVSS sources without a FIRST 
counterpart (~ 7.8 deg" 2 ; ~ 0.56 deg" 2 for sources 
brighter than lOmJy). For an astronomical ob- 
ject to be detectable by NVSS but not FIRST, it 
must be large and have a surface brightness too 
faint for the high-resolution FIRST survey On 
the other hand, it is possible that NVSS sources 
without FIRST counterparts are simply spurious, 
or have badly measured NVSS positions, i.e. their 
FIRST counterparts are outside the 30" matching 
radius. The top panel of Figure [5] shows the near- 
est neighbor distribution from the WENSS cat- 
alog for NVSS sources lacking a FIRST counter- 
part within 30". It peaks at small distances due to 
physically-associated objects, demonstrating that 
NVSS sources lacking a FIRST counterpart but 
with a WENSS match (within 55"), are dominated 
by real sources. If these are not spurious detec- 
tions, they belong to the category of extended, 
low-surface brightness objects missed by the high- 
resolution FIRST. 

Assuming that real NVSS sources have a coun- 
terpart in at least one of FIRST and WENSS, 
we estimate an upper limit for the fraction of 
spurious NVSS detections. Nearly 14% of NVSS 
objects in the region of survey overlap have no 
FIRST counterpart within 30", and 12.2% have no 
FIRST or WENSS counterpart, the latter within 
55". NVSS has a much fainter flux limit than 
WENSS (2.5 mJy as opposed to 18mJy): sources 
near the NVSS faint limit are likely to fall below 
the WENSS faint limit. For the sample of NVSS 
detections brighter than 18mJy, only 0.43% lack 
both FIRST and WENSS counterparts, which we 
interpret as an upper limit on the fraction of spuri- 
ous NVSS sources. The selection on which we base 
this estimate is biased against sources brighter at 
20 cm than at 92 cm; however, this type of radio 
source is rare, as we later show f *j4.2[) . Accord- 



ing to lde Vries et all (|2002l ). NVSS begins to lose 
completeness below 12mJy. We find that 2.3% 
of NVSS sources above this limit have no FIRST 
match. 
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3.6. Catalog availability 

The full catalog and several pre-made subsets 
are available for download on the catalog website^, 
which also provides a list and description of the 
data parameters. The datasets are described in 
Appendix [B] All 2,724,343 rows of the complete 
catalog (sample A) are available in a tarred archive 
of fits files, each covering a 5° wide strip in right 
ascension. A small subset of the catalog selected 
from « 100 deg 2 of sky, with 16,453 rows (sam- 
ple B), allows users to familiarize themselves with 
the data format. Several scientifically-useful sub- 
sets of the catalog are pre-made for convenience; 
some of these are analyzed in the next section. The 
first of these subsets contains radio fluxes and po- 
sitions for NVSS-FIRST associations (sample C), 
including 6 and 92 cm data, while another contains 
sources matched by all four radio surveys as well 
as the SDSS (sample F). Two more subsets con- 
sist of spectroscopic galaxies (sample G) and spec- 
troscopic quasars (sample H), detected by NVSS, 
FIRST, WENSS, and SDSS. Also provided are a 
sample of isolated^ FIRST-NVSS sources (sam- 
ple I) and isolated FIRST-NVSS-SDSS sources 
(sample J). Finally, a set of high-redshift galaxy 
candidates is available (sample K; see £|5.2[) . A 
more detailed description of these data files can 
be found in Appendix IbI 

4. PRELIMINARY CATALOG ANALY- 
SIS 

In this section, we investigate the distribution 
of sources in radio color-magnitude-morphology 
parameter space. We define the relevant mor- 
phology and spectral parameters using five ra- 
dio fluxes: peak FIRST flux (20 cm), integrated 
FIRST and NVSS fluxes (20cm), GB6 flux (6 cm), 
and WENSS flux (92 cm). While the unified cat- 
alog contains many more useful parameters inher- 
ited from the original catalogs, in this prelimi- 
nary analysis we focus only on these few. We use 
FIRST and NVSS flux measurements to define two 
morphology estimators, and use GB6 and WENSS 
fluxes with NVSS flux to compute radio spectral 
slopes. Morphology estimators are used to define 



three morphology classes, for which we construct 
radio "color-magnitude" diagrams. We also uti- 
lize optical photometric identification in the anal- 
ysis, classifying sources as SDSS point sources, ex- 
tended sources, or faint sources (non-detections). 

Hereafter our discussion covers samples selected 
using the conservative matching radii listed in Ta- 
ble [21 We settle upon these values by looking for 
an optimal balance between completeness and ef- 
ficiency of the matched samples, as illustrated in 
Fig. [3l The matching radii used herein are as fol- 
lows, with the larger values for all catalog asso- 
ciations given in parentheses: FIRST-NVSS, 25" 
(30"); FIRST- WENSS, 30" (120"); FIRST-GB6, 
70" (120"); FIRST-SDSS, 2" (60"). Estimated 
values of completeness and efficiency, calculated 
by comparing with random matching, are listed in 
Tabled 

Our analysis is limited to the 2955 deg region 
where the contributing surveys overlap (Fig. [I]) in 
order to control the selection criteria of the ana- 
lyzed samples. There are nearly 490,000 catalog 
entries in the overlap region, including multiple 
FIRST components matched to the same NVSS 
source, as well as duplicate entries (see H3.2p . To 
obtain a sample of physical FIRST-NVSS associ- 
ations, it is sufficient to select catalog entries con- 
taining an NVSS source and its nearest FIRST 
match: individual objects are typically not re- 
solved into multiple components in NVSS because 
of the survey's larger beam. This method selects 
~ 140, 000 unique, likely physical NVSS-FIRST 
associations. Of those, ~ 64, 000 are matched to 
a WENSS source, ~ 14, 000 to a GB6 source and 
- 48, 000 to an SDSS source. 

4.1. Morphology classification of radio 
sources 

4.1.1. "Simple" vs. "complex" FIRST-NVSS 
sources 

FIRST is a high-resolution interferometric sur- 
vey, and thus underestimates the flux of ex- 



tended and lobe-dom inated sources (jBecker et al 
Il995t iLu et~alll2007t) . Additionally, a multiple- 



3 http://www. astro. Washington, edu/akimball/radiocat/ 
4 Isolated here refers to an NVSS source with a single FIRST 
counterpart within 30". 



component source with core and lobes may be 
detected as three objects by FIRST but as only 
a single object in the lower-resolution NVSS. The 
difference between the two 20 cm magnitudes is 
thus a measurement of source morphology that 
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indicates angular extent and complexity. We de- 
fine 

At = £first — ^nvss- (2) 

Before applying At as a classifier we investigate its 
dependence on tNVSS for NVSS-FIRST pairs (sam- 
ple C) in Figure [6l The distribution in At is bi- 
modal, with the two peaks centered at At = and 
At = 0.7. Through the visual inspection of a thou- 
sand FIRST images (see below, <j4.1.3|) , we verified 
that sources in the At ~ locus are single compo- 
nent sources, while those in the At ~ .7 locus are 
multiple-component or extended. The lower panel 
shows the At distribution for four magnitude bins. 
We fit the sum of two Gaussians to each of the 
four distributions; best-fit parameters are listed 
in Table 03 The At distribution is similar for all 
magnitudes, although the overlap between the two 
peaks increases for fainter sources as flux measure- 
ment errors increase. We adopt At = 0.35 (solid 
horizontal line) as the separator between "simple" 
(At < 0.35) and "complex" (At > 0.35) sources. 

4.1.2. "Unresolved" vs. "resolved " FIRST 
sources 

The ratio of the two FIRST flux measurements — 
peak and integrated — yields a second measure of 
morphology: a dimcnsionlcss source concentration 
on ~ 5" scale. We define 

e=(f^) • (3) 

\ ± peak / 

Sources with 9 ~ are highly concentrated (unre- 
solved), while sources with larger 9 are extended 
(resolved). We adopt log(6» 2 ) = 0.05 (9 « 1.06) 
as the value separating resolved and unresolved 
sources. This choice is motivated by the distribu- 
tion of sources in the two-dimensional At vs. 9 
distribution, discussed below. 

4-1.3. Automatic morphology classification of 
radio sources 

The two-dimensional At vs. 9 distribution for 
sources with tNVSS < 13.5 (the brightest 30% 
of NVSS-FIRST pairs) is illustrated in Figure 
along with the marginal distributions for the two 
morphology parameters. The distribution in the 
lower panel suggests that 20 cm radio fluxes can 
be used to automatically separate FIRST-NVSS 



detections into three morphology classes: "com- 
plex" , (simple) "resolved" , and (simple) unre- 
solved or "compact" . 

This assertion was verified by extensive vi- 
sual inspection of FIRST images (Fig. [8]). Over 
1000 FIRST stamps (2x2 arcmin 2 ) of optically- 
identified radio quasars and radio galaxies were 
classified both visually and automatically (as 
above) into "complex", "resolved", and "com- 
pact" categories. A comparison of the results 
is given m Table [4] The two methods are con- 
sistent: over three quarters of all sources (81% 
of quasars and 76% of galaxies) receive identical 
classification from the two methods. Two initially 
surprising results are the significant fraction of ob- 
jects classified as "complex" by one method and 
"compact" by the other. In particular, 15% of 
visually-complex quasars were automatically clas- 
sified as "compact". As it turns out, the majority 
of these sources appear (by eye) to be asymmetric 
double-lobed sources with a very high flux ratio. 
Their FIRST images show two lobes, hence the 
visual "complex" classification. However, their 
total flux is due primarily to the brighter lobe, 
meaning FIRST and NVSS flux measurements 
are similar, resulting in an automatic classifica- 
tion of "compact" . Additionally, 26% of visually- 
compact galaxies were classified automatically as 
"complex" . An inspection of flux values reveals 
that many of these are borderline cases, with At 
value just above the At = 0.35 cutoff. It is likely 
that these galaxies have faint, diffuse emission 
which is detected by the NVSS but is not visible 
in the FIRST images. 

The automatic morphology classification presents 
difficulties because of non-unique pair matching, 
e.g. multiple FIRST detections matched to a sin- 
gle NVSS object. The resolved/unresolved clas- 
sification is based on FIRST fluxes and therefore 
applies to FIRST components individually. In 
choosing an analysis sample based on FIRST- 
NVSS pair matching, we retain only the closest 
FIRST match to an NVSS object (see Appendix[3 
for details) . The fluxes of the closest FIRST match 
are used, while more distant FIRST sources are ig- 
nored, when classifying a FIRST-NVSS pair as re- 
solved/unresolved. Approximately 17% of NVSS 
detections in the FIRST footprint have multi- 
ple FIRST matches within 30". Most of these 
sources (88%) are classified as complex, presum- 
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ably because the NVSS flux encompasses multiple 
FIRST components. We therefore apply the re- 
solved/unresolved classification only to "simple" 
sources, resulting in the three morphology clas- 
sifications. (In other words, we do not separate 
the complex sources into "complex resolved" and 
"complex unresolved".) 

Figure [9] shows the magnitude distribution 
for NVSS-FIRST associations of each morphol- 
ogy class, and the fraction in each class of all 
sources brighter than a limiting £nvss magnitude. 
There are similar numbers of complex and com- 
pact sources at all magnitudes; there are fewer re- 
solved sources, and the fraction increases slightly 
at fainter flux values. Overall, the fractions re- 
main roughly constant with flux. In the next 
section, we investigate correlations between radio 
morphology and radio spectral slope. 

4.2. Radio color-magnitude and color- 
color diagrams 

In the previous section, we described an auto- 
matic morphology classification for NVSS-FIRST 
associations based on their 20 cm fluxes. The ad- 
dition of WENSS (92 cm) and GB6 (6 cm) data 
allows us to compute radio spectral slope, and to 
study its correlation with radio morphology and 
magnitude. 

For two magnitudes t\ and £2, we compute a 
spectral index a, defined by F v oc v a , as 



0.4 



H 2 



log(A 4l /A t2 ) 



(*i -*a). 



(4) 



The value a describes the average spectral slope 
between wavelengths and Xt 2 ; it is proportional 
to the color t\ —t^. Therefore 



4° = 0.765(i N vss 



tGB6), 



a| 2 = 0.337(i\vENSs — *gb6)j 



(5) 
(6) 



and 



a 9 2 2 = 0.604(i WENSS - iNvss ) ■ (7) 

Note that we use NVSS rather than FIRST when 
calculating spectral indices for the analysis pre- 
sented here, as NVSS is more likely to detect all 
the 20 cm flux from a source due to its larger beam- 
size. 

Typical extragalactic radio sources have spec- 
tral index in the range — 1 < a < in the 



radio regime (e.g., ICondonl Il984l ; iBlundell et al 
1999), with a = —0.5 as a typical value separating 



steep-spectrum and fl at-spectrum sources (e.g., 
Urrv fc Padovanil ll995). Emission from extended 



radio lobes, such as those seen in classical FR2 ob- 



jects ([Fanaroff k. RilevHl974l ) tends to be steeper, 
with the steepening attributed to synchrotron, in- 
verse Compton, and adiabatic losses acting on the 
relativistic electrons. The flatter emission from 
compact quasar cores /jets is in terpreted as self- 
absorbed sy nchrotron emission (lKroliklll999l) . As 
described in iJarvis fc McLurd (|2006h . radio spec- 
tral index is thought to correlate with source orien- 
tation in AGN unification schemes: sources viewed 
along the jet axis show flat radio spectra, sources 
viewed along the dust torus show steep radio spec- 
tra (since only the lobes are visible), and inter- 
mediate orientations result in intermediate spec- 
tral slopes due to mixing of the flat- and steep- 
spectrum emission. A minority of observed radio 
sources fall into the gigahertz-peaked-spectrum 
(GPS; 10%) and compact-stee p-spectrum (CSS; 
30%) categories (lO'Deal Il998h . GPS and CSS 
sources are similar to smaller versions of FR2s, 
but have convex spectra that peak around 1 GHz 
in the observer framed 

Two significant systematic issues could cause 
errors in our measure of radio spectral index. The 
first of these is the possible variability in radio 
flux over the timescale of the survey observations: 
GB6 observations were carried out in 1986 and 
1987, WENSS operated in the early 1990s, and 
NVSS observations were taken in the first half 
of 1998. AGN have shown flux variations over 
all radio wavelengths on timescales ranging from 
less than a day to years. Intrinsic variability is 
thought to be due to jet v ariations, such as shocks 
propagating along t he jet ( Marscher fc Gear! 19851: 
iHughes et al.l Il989h , with Doppler boosting act- 
ing to exaggerate the effect. Therefore, variabil- 
ity should be more prevalent in sources viewed 
along or near the jet axis, which is consistent 
with observations: blazars, quasars, flat-spectrum 
sources, and compact sources show higher vari- 
ability among a greater fraction of sources than 



5 GPS galaxies are thought to be young, evolving radio 
sources that expand through a CSS stage and eventually 
turn into FR1 or FR2 galaxies llO'Dealll99l l de Vries et al.l 
2007). GPS and CSS sources that peak at sufficiently short 
wavelengths can have positive a in the unified catalog. 
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contro l samples. For example, iBarvainis et al 
([20051 ) observed variability at the 5-8% level over a 
two-year timescale for core-dominated quasars ob- 
served at 8.4 GHz. They also observed some large 
variations on month-long timescales, but only one 
in five of their sources varied at the 10-30% level, 
and most vari e d by less than 20%. Meanwhile, 
de Vries et ail <|2004h found that roughly 2% of 
their sample of unresolved FIRST sources show 
significant variability (> 4cr). In a sample of 
3600 sources brighter t han 0.1 mJy at 1.4 GHz, 
Rvs fc Machalskil (|l990h found that less than 1% 
were variable, and that all of those have flat radio 
spectra, with 86% being core-dominated sources. 
As far as measurements for individual sources are 
concerned, a 10% flux variation at one of the two 
relevant wavelengths corresponds to a change of 
0.08 in spectral index, while a 30% flux varia- 
tion corresponds to a change of about 0.2 in spec- 
tral index. Results from the literature discussed 
above suggest that about 1% of the unified cat- 
alog sources varied significantly over the range of 
survey observations, and that the variation tended 
to be at the 10% level. Those sources are likely to 
be concentrated in the compact radio morphology 
sample with flat radio spectra (a > —0.5). This 
variability contributes to scatter in the distribu- 
tion of a. 

The second issue which may cause erroneous 
calculations of spectral index is the fact that the 
radio surveys have different spatial resolutions. 
WENSS and NVSS both have ~ 50" resolution 
(see Table [TJ, while GB6 has a resolution of about 
200", and therefore gathers flux from a much 
larger area. This issue will effect spectral in- 
dex for radio sources resolved into separate NVSS 
components, such as a double-lobed radio galaxy 
with lobe separation ~ 100". If a single NVSS 
component contains only half of a source's total 
20 cm flux while GB6 sees all of the 6 cm flux, 
then a§° will be too shallow by almost 0.6 for that 
source. We note that since WENSS and NVSS 
have similar angular resolution, this issue is un- 
likely to significantly affect the o^o measurements. 

Figure [10] presents color-magnitude and color- 
color diagrams for the three radio morphology 
classes detected in all four radio surveys (listed as 
sample E in Appendix iBj) . Similar color-color di- 
agrams to th ose in the m i ddle co lumn were first 
presented in llvezic et"aL (|2004al Fig. 2), with 



slightly different selection criteria. The middle col- 
umn demonstrates that the majority of sources in 
all three morphology classes have spectral indices 
a 20 ~ a l° K ~0-8. This resu lt is in agreement 



with earlier observations (e.g., Kellermannl 11964 ; 



Laing et al1ll983t Zhang et al.ll2003 ). and suggests 
that most strong radio sources have a fairly con- 
stant power-law slope from 6 to 92 cm. However, 
the radio compact sample (lower middle panel) 
has a significant fraction of flat-spectrum sources 
(a ~ 0) as well. In the next section, we investigate 
the two groups of compact sources using optical 
identifications, and show that they are due to a 
population of galaxies and a population of optical 
point sources, with different distributions in radio 
color-color space. Note that in the color-color dia- 
grams, the WENSS faint limit biases against large 
positive a, while the GB6 limit biases against large 
negative a. The strong correlation between 
and (£nvss) in the right column is due to the GB6 
survey sensitivity: only very bright steep sources 
are detected by GB6. A small fraction of sources 
(< 2% with £nvss = 12) have a negative value for 
a|o but a positive value for a|°, i.e. are fainter 
at 20 cm than at 92 cm and 6 cm. Such sources 
either have highly unusual radio spectra, or varied 
significantly in brightness over the range of obser- 
vations. 

We quantify the fraction of compact sources 
with steep and flat spectra in Figure [TTJ The 
sources are divided into four subclasses depending 
on their spectral shape between 6 cm and 92 cm. 
The top panel shows an example shape for each 
subclass. The "steep" and "flat" designations re- 
fer to sources with monotonic radio flux measure- 
ments. The middle panel shows the magnitude 
distribution of each subclass for compact sources 
with £nvss < 12. Slightly less than two thirds of 
the compact sources are steep-spectrum, whereas 
94% of complex and resolved sources are steep- 
spectrum. We find that one quarter of compact 
sources are fiat-spectrum or peaked between 6 cm 
and 92 cm; these values are consistent w ith the 
fractions of GPS and CSS sources found bv lO'Deal 
fjl998h . The fractions of steep, flat, and peaked 
sources remain nearly con stant with magnitude 
(bottom panel). However, Ide Vries et al.l ([2002) 



observed that the ratios change at fainter fluxes, 
where the fraction of flat-spectrum sources rises 
sharply (their Fig. 1). Note that because of the 
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flux limits of the contributing radio surveys, we are 
more sensitive to inverted-spectrum sources than 
to peaked-spectrum sources. 

4.3. Optical identification of radio sources 

The addition of SDSS data enables the sep- 
aration of radio sources into galaxies (i.e., ex- 
tended optical sources), optical point sources, and 
faint optical sources (i.e. undetected by SDSS). 
Combining the three optical categories with the 
three radio morphology categories yields nine opti- 
cal/radio classes. As discussed in §2.2. It the sam- 
ple of optical point sources may include a small 
number of unresolved galaxies and radio stars in 
addition to radio quasars. 

4.3.1. The NVSS-FIRST- WENSS- GB6 sample 

We begin by discussing the approximately 
12,000 sources detected by all four radio surveys 
(required in order to measure the two radio spec- 
tral indices a^o and a|°). Over 38% of these have 
an SDSS counterpart within 2". This is slightly 
larger than the fraction of FIRST sources with 
an optical counterpart (32%), and may be due to 
a bias towards the bright end induced by requir- 
ing GB6 and WENSS detections. This dataset is 
listed as Sample E in Appendix [BJ the optically- 
matched subset is Sample F. 

Figure [12] shows £nvss magnitude distributions 
for sources detected by all four radio surveys and 
brighter than tNvss = 12, corresponding roughly 
to the GB6 detection limit (see §3.5]) . For the 
optically-matched subset, most sources have com- 
pact radio morphology, while the fraction of com- 
plex sources drops dramatically. This is partly 
due to the fact that radio-optical matche d samples 
are b iased against double- lobed sources ( Lu et all 
2007|). The fractions of sources in each morphol- 
ogy class remain roughly constant with iNVSS even 
for the optically-matched subset. In all three radio 
morphology classes, galaxies and point sources are 
found in roughly equal ratios at all magnitudes. 

Figure [12] shows radio color-color diagrams for 
the nine radio morphology/optical identification 
classes. This figure is an expansion of the middle 
column of FigurefTOl with the addition of the three 
optical identification categories. Faint sources, 
galaxies, and point sources show different color- 
color distributions in the three morphology classes, 



particularly in the compact sample. For the com- 
plex sample, faint optical sources peak at steep 
values of a ~ — 1 while point sources peak near 
a ~ —0.5. Galaxies are more tightly clustered 
around the peak than point sources. The sharpest 
separation in behavior is seen in the compact mor- 
phology sample (bottom row). The optically- faint 
and galaxy samples tend to have steep spectral 
indices (a < —0.5), while point sources tend to 
have flat spectral indices (a > —0.5). This con- 
firms, with much improved statistics, the long- 
known distinction between flat spectrum radio 
quasars which are dominated by core emission and 
steep spectrum radi o quasars wh ich are dominated 
by lobe emission (Krolik 1999). The compact 
point sources are also bimodal, with both a flat- 
spectrum and a steep-spectrum peak, while point 
sources in the complex and resolved morphology 
classes are primarily steep-spectrum. The steep- 
spectrum peak in the compact morphology class 
(lower right panel) may be due to distant, un- 
resolved galaxies contaminating the point source 
sample. 

The separation of the steep- and flat-spectrum 
peaks in color-color space is maximized along the 
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- a?g locus. We thus take the average of 
the two indices to quantitatively investigate spec- 
tral slope distributions, as presented in Figure [14] 
Complex and radio resolved classes contain pri- 
marily steep spectrum sources, while the optically- 
identified radio compact sample is bimodal with a 
peak due to flat-spectrum point sources and a peak 
due to steep spectrum galaxies. The optically- 
faint subset of the compact morphology class is 
biased toward steep-spectrum sources (compare 
panels A and C). The distributions for the op- 
tically faint sample are similar to those of the 
galaxy sample. These histograms are based on the 
largest multi-wavelength samples of radio sources 
ever constructed. Table[5]compares the mean, me- 
dian, and FWHM for the distributions in panels 
A, D, E, and F. 

4.3.2. The NVSS-FIRST- WENSS sample 

The majority of radio sources have a < —0.5, 
i.e. they are brighter at longer radio wavelengths. 
The sample of sources detected by all four radio 
surveys is therefore flux limited by the short wave- 
length GB6 survey. Dropping the requirement of 
a detection in GB6 increases the number of se- 
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lected sources to 64,000, a factor of five, while 
spectral slope information is still available from 
a|o ■ The a^o distributions for the larger sample of 
NVSS-FIRST-WENSS detections (not shown) in- 
clude much fainter sources; the observed bimodal- 
ity seen in Figure [M] is still present, although the 
separation is weaker. Mean, median, and FWHM 
of the alo distributions for the samples in Fig- 
ure [H are shown in Table [6] We use the NVSS- 
FIRST-WENSS subsample (listed as sample D in 
Appendix [B| for all further data analysis. 

4-3.3. Summary of source characteristics in the 
unified radio catalog 

Figure [15] shows radio color-magnitude dia- 
grams for NVSS-FIRST-WENSS sources from 
the catalog overlap region, providing a qualita- 
tive summary of fundamental properties of radio 
sources in the unified catalog. This figure ex- 
pands the left panel of Figure [TOl by adding classes 
defined by the optical identification categories. 
However, Figure [10] was limited to sources ob- 
served by GB6 in addition to NVSS, FIRST, and 
WENSS, resulting in a sample of only the bright- 
est radio sources. Figure [15] demonstrates that 
GB6 detections are not necessary to obtain an 
extensive characterization of these radio sources, 
including radio spectral slope. Sources as faint 
as £nvss ~ 13-9 can be studied without biasing 
toward steep-spectrum sources. The sample pre- 
sented here contains 63,660 radio sources (optical 
subset: 18,728). Constraining to £nvss < 13-9 re- 
sults in a sample of 50,046 (optical subset: 15,424). 
This catalog represents a significant increase in 
comprehensive radio data available to the scien- 
tific community. 

Figure [15] provides a visualization of the uni- 
fied radio catalog in a five dimensional param- 
eter space defined by four radio magnitudes 
(*NVSS,*first (integrated), and t F i RST (peak) at 
20 cm; £wenss a t 92 cm) and one optical morphol- 
ogy parameter (optical light concentration; see 
i j2.2.1|) . The four radio magnitudes are combined 
into parameters indicating magnitude, morphol- 
ogy, and spectral slope. Using three radio mor- 
phology classes, radio magnitude, spectral slope, 
and three optical identification/morphology sub- 
classes, we have summarized the five-dimensional 
parameter space using the nine two-dimensional 
panels (radio color-magnitude diagrams) in this 



figure. 

The marginal distributions of spectral index 
«2o f° r t ne sample in Figure [15] are similar to the 
distributions shown in Figure RBI However, the bi- 
modality seen in the compact sample is reduced. 
This is due both to the inclusion of fainter sources, 
and the fact that one must use a^o? rather than 
the average of q^q and a§°, to describe the radio 
spectra of a sample that does not require detection 
at 6 cm. 

Magnitude distributions for the nine radio/optical 
morphology subclasses are shown in Figure 1161 
This is the equivalent to Figure [12] but corre- 
sponding to the five times larger sample of NVSS- 
FIRST-WENSS sources. This figure shows the 
striking result that, at magnitudes fainter than 
fervss = 11-5 (wlOOmJy), the density of opti- 
cal point sources flattens out dramatically, rather 
than continuing to rise as would be expected for 
a homogenous source distribution. This catalog 
is the first sample of radio sources both large 
and deep enough to investigate such faint source 
counts. The same behavior is seen in all morphol- 
ogy classes: point source counts flatten out while 
galaxy counts continue to rise at faint magnitudes. 

4.4. Optical analysis aided by radio infor- 
mation 

In the previous section, we showed that analy- 
sis in radio space is greatly aided by the addition 
of optical information. Similarly, the standard op- 
tical analysis, such as determination of optical lu- 
minosity function and its evolution with redshift, 
can benefit from the addition of radio observa- 
tions. For example, one can ask whether radio 
loudness and radio spectral index correlate with 
optical quantities (e.g., "Can SDSS data alone pre- 
dict the strength and spectral properties of radio 
emission?"). In this section, we take initial steps 
in this direction by examining the dependence of 
radio spectral slope on redshift and optical mag- 
nitude for quasars, and on optical magnitude and 
color for galaxies. We also investigate the relation- 
ship between radio loudness, radio spectral slope, 
and optical color. 

We adopt the spectral index af 2 as a measure 
of radio loudness, where 

af = 0.0657(%enss-«)- (8) 
This spectral index relates the 92 cm flux to the 
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SDSS i band flufi and can be converted to the 
often used r adio loudness mea sure, i?,: = 0A(i — 
iFiRST) (e.g. Ilvezic et al.ll2002h via 



a\ z = O.llo^ - 0.16 B4 



(9) 



A minimum in the ra dio loudness distribu tion at 
Ri ~ 1 reported by Ilvezic et all (|2004bl ) corre- 
sponds to 
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O.lla^- 0.16 ~ -0.25. 



(10) 



The NVSS-FIRST-WENSS sample analyzed here 
is not sufficiently sensitive to investigate radio 
loudness bimodality, due to the WENSS flux limit. 

Throughout this section, we focus on sources 
observed by the NVSS, FIRST, WENSS, and 
SDSS surveys, with optical spectra from the 
SDSS. We assume a cosmology with Hq = 
70 km s^Mpc" 1 , M = 0.3, and fl A = 0.7. 

4--4-1- Spectroscopically-identified quasars 

To find a sample of spectroscopically-identified 
radio quasars, we match unified catalog sources to 
the F ifth Data Release (DR 5) SDSS Quasar Cat- 
alog (jSchneider et al.ll2007t hereafter S07). This 
catalog is a cleaner quasar sample than sources 
identified as quasars by the SDSS spectroscopic 
software pipeline: the S07 spectra have been visu- 
ally examined to verify that they have at least one 
broad emission line (FWHM > 1000 km s _1 ) or 
are unambiguously broadline quasars. Candidates 
for the visual examination included all sources se- 
lected by the SDSS quasar target selection algo- 
rithm, and all sources identified as quasars after 
spectroscopic observation. 

The sample used here consists of the 1,288 
NVSS-FIRST-WENSS sources matched to an S07 
source within 2" (sample H). Note that because 
the quasars were discovered via various SDSS tar- 
get selection algorithms, they do not represent a 
complete statistical dataset. Indeed, proximity 
to a FIRST source is one of the selection meth- 
ods. One therefore must be careful when creating 
statistical samples from this data set; e.g. one 
must choose quasars identified by the same tar- 
get method, and account for sources in regions of 



3 An alternative m easure of rad i o loud ness is based purely on 
radio luminosity. Ilvezic et al] 1 20021 found that the two dif- 
ferent quantities yield similar radio loudness classifications, 
owing to strong selection effects in flux-limited samples. 



the sky where different target selections were used. 
Note however that the SDSS optical quasar sam- 
ple is estimated to be complete at the 94.6% level 
for the magni tude range 16 < i < 1 9. and redshifts 



less than 5.8 (Ri chards et~aTll2002l l. 



Figure [17] investigates the dependence of ra- 
dio spectral slope on optical luminosity and red- 
shift for quasars with different radio morphologies. 
Optical K-corrections used to calculate absolute 
magnitude Mj assume an o ptical spectral index 
a = —0.5 (jlvezic et al]|2002l) . Panel A shows that 
the radio spectral slope for radio quasars is inde- 
pendent of redshift and is a strong function of ab- 
solute magnitude: quasars with the steepest radio 
spectra tend to be optically faint. We investigate 
this correlation as a function of radio morphol- 
ogy in panels B-D: the trend is seen quite strongly 
in the compact morphology class, but is weak for 
the complex and resolved sources. Additionally, 
compact quasars have the flattest radio spectra, 
while radio spectra of complex sources are rela- 
tively steep. 

A more quantitative visualization of the cor- 
relation between radio spectral slope and optical 
luminosity is shown in Panels E and F, for ra- 
dio complex and radio compact sources, respec- 
tively. The median value in unit magnitude bins 
is indicated by black circles. The median spec- 
tral slope for complex sources is constant with 
luminosity, while the median slope for compact 
sources is a strong function of luminosity. This 
trend may indicate a physical change in radio 
spectral slope with optical luminosity, an orien- 
tation effect, and/or intrinsic variation in spec- 
tral slope combined with luminous sources be- 
ing observed at higher redshifts (i.e., a radio K- 
correction effect). Figure [18] shows as an ex- 
ample a simulated spectrum that steepens for 
A > 20 cm, which would show the observed trend 
when viewed at different redshifts. It has been 
suggested however that spectral slope is an in - 
dication of orientation ( Jarvis fc McLurd[2006l) . 
and therefore of increased Doppler beaming. The 
spectrum of a quasar viewed along a line of sight 
close to the jet axis is highly core-dominated. It 
would thus appear flatter than a source viewed 
off-axis, while beaming would lead to a higher 
observed luminosity. Complex sources (Panel E) 
can be expected not to show the same trend, 
as they almost certainly have extended emission 
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from steep-spectrum radio lobes, biasing against 
sources with strongly-beamed jets. On the other 
hand, the compact sample may contain unresolved 
core-lobe sources, especially at higher redshift. 
An angular size of 5" (FIRST resolution) cor- 
re sponds to ~ 40kpc at reds hift 1. According 
to iBhmdell k Rawlingsl (|2000l Fig. 2), a typical 
quasar can reach this size in ~ 10 6 — 10 7 years 
while the maximum active lifetime is estimated to 
be 10 8 — 10 9 years. However, a line-of-sight align- 
ment near the jet axis would significantly shorten 
the projected linear size, decreasing the likelihood 
of resolving a core- lobe source. Note that the com- 
bination of a flat core component and a steep lobe 
component can result in a spectrum of the type 
shown in Figure[18] The trend observed in Panel F 
may therefore be due to an intertwined, and insep- 
arable, combination of the above effects. 

Equivalent plots to Panels E and F (not shown) 
indicate that the median value of a§° at all opti- 
cal luminosities corresponds to a flat spectrum be- 
tween 6 and 20 cm in all three morphology classes, 
which lends support to the K-correction argu- 
ment. There are quasars whose spectra are con- 
sistent with the simulated example shown in Fig- 
ure [TS1 but n ot all , as demonstrated in Fig. 2 of 



Kuehr et alj (|198lh . However, the resolution of 
their observations was very low. To accurately de- 
termine the radio spectral shape of quasars will 
require better frequency resolution, which will be 
possible with the future Expanded VLA project 
(lNapieil[2006h . 

The dependence of radio loudness on radio 
spectral slope and optical color for quasars is in- 
vestigated in Figure 1191 Sources are divided into 
steep- (left) and flat-spectrum (right). We re- 
quire twENSS < 12.9 so as not to bias against any 
sources with a|o ~ —1-5. We require i < 19 to 
ensure a complete sample of quasars. Based on 
those two limits, we have selected a region of ra- 
dio/optical magnitude space which does not bias 
against radio quiet or radio loud sources (top two 
panels). The radio loudness distribution of steep- 
spectrum quasars (lower left panel) is strongly de- 
pendent on color: a Kolmogorov-Smirnov test re- 
jects the null hypothesis that the two histograms 
share the same parent distribution at the 99.6% 
level. No such dependence is seen for the flat- 
spectrum quasars (lower right panel). This paper 
presents the first observation of this effect. Note 



that the requirement of a detection at 92 cm biases 
these samples to bright radio fluxes, and prevents 
us from reaching the expected local minimum in 



the oa distribution at oa 



-0.25. 



4-4-%- Spectroscopically-identified galaxies 

The sample of spectroscopically-identified galax- 
ies (sample G) is the set of 2,885 NVSS-FIRST- 
WENSS-SDSS matches whose spectra identify 
them as galaxies, according to the SDSS auto- 
matic classification software. Here, we restrict our 
analysis to a flux-limited sample with r < 17.77, 
in order to ch oose galaxies from t he SDSS main 
galaxy sample ( Strauss et al. 20021 ). The resulting 
flux-limited sample consists of 1,656 spectroscopic 
galaxies. 

Figure [201 investigates the dependence of radio 
spectral slope on optical luminosity and color for 
galaxies. Complex galaxies tend to have steeper 
radio slopes than compact galaxies. Radio sources 
with flat spectrum emission are generally under- 
stood to be qua sar cores (i.e ., quasar sources with- 
out radio lobes, lKroliklll999h . Quasar core sources 
are compact in the radio due to their small size and 
typically large distances. Sources in panel D there- 
fore beg the question: what are the flat-spectrum, 
radio compact sources that are optically-identified 
as galaxies? It is possible that those sources are 
quasars too weak in the optical to overpower emis- 
sion from the galaxy, observed as galaxies with 
AGN (e.g. Seyfert galaxies, LINERS), a possibil- 
ity we will investigate in future work. 

Figure [2H explores the dependence of radio 
loudness of galaxies on optical luminosity and 
color. The top row shows the distribution of steep- 
(left) and fiat-spectrum (right) sources in optical 
color-magnitude diagrams with points color-coded 
by a 92 . Note that only very luminous galaxies are 
radio loud (a? 2 < -0.35 for M r < -22). The 
lower two panels show the distribution of a 9 f for 
steep and flat radio galaxies divided by u — r color. 
Similar to Fig. 1191 we focus on an unbiased re- 
gion of radio/optical magnitude space, based on 
the limits iwENSS < 12.9 and i < 17. The lower 
left panel shows a significant difference in radio 
loudness for blue (dashed line) and red (solid line) 
steep-spectrum radio galaxies. The blue galaxies 



have a median a? 2 of 
dian value for red galaxies is 
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-0.13, while the me- 
-0.26. A dif- 



ference of 0.13 in a 92 (Eq. [8]) corresponds to a 
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change of about 2 magnitudes, i.e. steep-spectrum 
red galaxies are about two orders of magnitude 
brighter in the radio than steep-spectrum blue 
galaxies of the same optical brightness. The dis- 
tributions of radio loudness for blue and red flat- 
spectrum galaxies are not significantly different. 

5. CATALOG APPLICATIONS/FUTURE 
WORK 

The information provided by this catalog will 
enable many diverse studies. In this section, we 
briefly remark on applications toward the investi- 
gation of the radio quasar/galaxy unification the- 
ory, the search for radio stars, and the selection of 
possible high-redshift galaxies. We then describe 
the usefulness of the catalog for studies of the evo- 
lution of the radio universe by comparison with 
models, which we discuss in detail in a companion 
paper (A. Kimball et al. in preparation). 



objects will lend evidence for or against the uni- 
fication paradigm by the investigation of num- 
ber statistics, the size versus orientation distribu- 
tion, and environment. The unification scenario 
has successfully explained many observed proper- 
ties of bright radio sources in previous research 
dBarthe]|ll989HUrrv et alJll99ll [Padovani fc Urrv 



19921 iLister et al.lll994l lUrrv fc Padovanil ll995j T 

but recent studies suggest there may be some 
intrinsic differences betwee n radio quasars and 
galaxies (|Willott et al.ll2002f) . 



5.1. 



Unification paradigm for radio-loud 
active galactic nuclei 



The orientation-based unification scenario (Urry_& 



19951 ) was originally motivated by the similarity 
of radio emission from sources that appear very 
different optically, i.e. quasars and galaxies. The 
theory assumes that sources whose emission is 
dominated by a radio core or lobes are members 
of the same "parent" population, but differ in ap- 
pearance because their highly anisotropic emission 
is viewed from different angles. This anisotropy 
is predominantly due to relativistic motion of the 
plasma in the inner jets and Doppler boosting 
("beaming") when the angle between the line 
of sight and the plasma velocity vector is small 
(lKrolikl [r999). Hence, the same object could ap- 
pear as a core source when viewed at small angles 
(with the "boosted" core outshining the lobes), 
and as an extended double-lobe or a core-lobe 
source otherwise. Understanding the radiation 
anisotropics in AGNs is required to unify the 
different types; that is, to identify each single, 
underlying AGN type that gives rise to differ- 
ent observed classes through different ori entations 
(|Antonuccilll993HUrrv fc Padovanilll995l) . 

The unified catalog allows users to classify 
many quasars and galaxies according to their ra- 
dio morphology (e.g., core sources, lobe sources) 
and spectral behavior. Statistical studies of such 



5.2. A method for selecting high-redshift 
(z > 1) galaxy candidates 

Galaxies at high redshift are important for 
studies of large scale structure and galaxy evolu- 
tion, and in fact the higher the redshift of a galaxy 
the more useful it is for these studies. Candidate 
high-redshift galaxies can be selected from the uni- 
fied catalog based on their radio morphology, ra- 
dio spectral slope, and lack of an optical counter- 
part. The lack of optical counterpart requirement 
is int ended to select sources so distant that their 
-ESid^c^Jii-vcd optical emission lies blueward of the rest- 



frame 4000A break (e.g., iMadau et ai1ll996l ). As 
the optical sources used herein were selected in 
the r and i bands at 6 165 A and 7481 A respec- 
tively, this requirement will tend to select galax- 
ies at redshifts of z ~ 1 and higher. A steep ra- 
dio spectrum has often been used as a criterion 
to find high-redshift sources, owing to the fact 
that higher redshift objects are observed at higher 
rest-frame frequencies, and because radio galaxy 
spectra tend to flatten below ~ 3 00 MHz in their 
rest frame (e.g., Cruz et al. 20061 and references 
therein). Additionally, as discussed in ^14.31 in the 
compact morphology subclass a steep radio spec- 
trum is a likely indicator of an optically-resolved 
source. 

To find a sample of candidate high-redshift 
galaxies, we select sources which are unresolved by 
FIRST, undetected by the SDSS, and have steep 
radio spectra. Specifically, we begin with the sam- 
ple of objects identified by NVSS, FIRST, and 
WENSS which have compact radio morphology. 
We require q^q < — 0.5 and no SDSS match within 
3". Note that heavily dust-obscured galaxies are 
also likely to be found in such a sample. In the cat- 
alog overlap region, the above criteria select 9,953 
sources. Visual examination of SDSS images of all 
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the candidates revealed 334 suspect sources, i.e. 
radio positions coinciding with a diffraction spike, 
scattered light, a nearby star or galaxy, etc. Re- 
moving the suspicious sources results in a sample 
of 9,619 candidate z > 1 galaxies. It will be inter- 
esting to compare this sample with deeper optical 
surveys and observations at other wavelengths, as 
the data become available. This sample is avail- 
able for download on the catalog website (sam- 
ple K). 

5.3. The search for radio stars 

With the great increase in the sensitivity of ra- 
dio surveys in the last several decades, along with 
the more accurate source positions afforded by ra- 
dio interferometry, comes the ability to search for 
fainter and rarer radio sources. A small fraction 
of stars have significant non-t hermal radio emis- 
sion, such as dMe flare stars ( White et al J 119891 
and cataclysmic va riables ( Chanmugaml Il987t 
Mason fe Gravl 120071 ) , but most stars have only 
weak thermal radio emission. Because quasars 
and stars lie in different locations in SDSS color- 



color diagrams (jlvezic et al . 2002; Richards et al 



20011 ). it may be possible to select a clean sample 
of radio stars based on their photometric colors. 
To investigate this possibility, we select a sample 
of radio star candidates from the unified radio 
catalog. 

We begin by selecting point sources from the 
catalog overlap region which coincide with a 
FIRST radio detection. We use a conservative 
matching distance of 1" to limit the contamina- 
tion due to random matches. To ensure a sam- 
ple for which the majority have spectra, we re- 
quire i < 190 We eliminate the region of color- 
color space where quasars are commonly found: 
u-g < 0.8 and -0.2 < g-r < 0.6. 

The number of candidates remaining after each 
step of the selection algorithm is shown in Ta- 
ble [7l The resulting sample of 532 candidate ra- 
dio stars contains 406 objects with SDSS spec- 



7 The SDSS quasar target selection algorithm flags all 
% < 19.1 sources within 2" of a FIRST detection 
l lSchneider et al.|[2007b . 



trail From the SDSS spectral classifications^ 
we determine that 78 are stars, 313 are quasars, 
and 15 are galaxies. The results of applying 
the selection to the SPSS DR5 quasar catalog 
(jSchneider et al.ll2007l ) are also shown for compar- 
ison. T he completeness of the quasar catalog is 



95% ( Richards et all 12002; Vandcn Ber k et al 



2005) . Seven of the known quasars were identified 



as galaxies by the SDSS spectroscopic pipeline, 
and two as stars. The latter are actually the super- 
position of a star and a quasar: the spectra are pri- 
marily stellar, but also contain at least one broad 
emission line. Note that 33 of the candidate ra- 
dio sources classified as "quasar" are absent from 
the quasar catalog. These are sources that were 
spectrocopically-observed for the sixth SDSS data 
release, but are not in DR5. 

The sky density of unresolved i < 19 SDSS 
sources is ~ 1000 deg~ 2 , while the density of 
quasars is ~ lOdeg -2 . However, out of nearly 6 
million point sources with i < 19, we found only 
78 confirmed radio stars, implying that the like- 
lihood of a star being a radio source is approxi- 
mately 10~ 5 , while the likelihood for a quasar is 
about 8 x 10~ 2 . Consequently, despite using a con- 
servative color cut to eliminate probable quasars, 
the sample is dominated by radio quasars rather 
than radio stars. Using this color-color selection is 
therefore not an efficient way to find radio stars. 
One must have spectral classification in order to 
select a clean radio star sample; photometric ob- 
servations alone are not sufficient. 

5.4. Radio galaxy evolution models 

In the companion paper, we attempt to un- 
derstand and quantitatively explain the observed 
source counts for the nine subclasses defined 
by radio morphology and optica l identification 
(Fig. [13|) . with the aid of models. iBarai fe Wiita 



8 The remaining 126 stellar candidates do not have spectra 
for one or more of the following reasons: they are not in 
the region of sky included in the SDSS DR6 spectroscopic 
sample, they are too close on the sky to another spectral 
target, they failed the quasar target selection algorithm 
because of the bright limit or bad photometry, or they were 
not targeted because the photometry in a previous SDSS 
data release led to an extended rather than point source 
optical morphology identification. 

9 Spectroscopic identifications determined using the SDSS 
"SpecBS" Princeton reductions. For details, see 
http: / /spectro. princeton.edu/ 1 
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( 20061 2007 ) investigated a large number of in- 
put parameter variations for several radio galaxy 
evolution models and compared the resulting dis- 
tributions of radio power, spectral index, linear 
size, and redshift t o the low- frequency f~2m] 
Cambridge catalogs (ILaing et al ll983tlEaleJl98a 



McGilchrist et al.Nl990h . Thev found that none of 
the models was able to provide good simultaneous 
fits to all observables for all three catalogs, but 
noted that the observational datasets were too 
small to adequately constrain the model parame- 
ters. 

The catalog presented in this paper represents 
a significant increase in the size of radio samples 
for which luminosity, size, spectral slope, and red- 
shift can be measured. We will compare the ob- 
served distributions from the unified catalog to 
distributions of sources from mock catalogs (ob- 
tained by applying observational selection crite- 
ria to the computer-generated radio skies). The 
distributions will help to constrain the increas- 
ingly sophisticated models of radio galaxy evolu- 
tion. The models will thus help explain the physics 
and evolution leading to the distributions seen in 
the unified catalog. 



5.4.1. 



Defining 
sample 



a long-wavelength flux-limited 



The radio galaxy evolution models currently 
perform best at long wavelengths, where emission 
is dominated by the radio lobes. The most ro- 
bust approach for comparing models and obser- 
vations is therefore to define a flux-limited sam- 
ple using long-wavelength measurements, such as 
the 92 cm flux from the WENSS survey. Fortu- 
nately, because most radio sources are brighter 
at longer radio wavelengths, a source observed in 
WENSS (faint limit of 18mJy) is almost guaran- 
teed to be detected at 20 cm (faint limits of 1 mJy 
for FIRST and 2.5 mJy for NVSS). However, the 
92 cm flux limit must be chosen carefully in order 
to avoid biasing against sources with radio spec- 
tra so steep that they are too faint at 20 cm to 
be detected by NVSS. The color-magnitude dia- 
grams and spectral slope histograms presented in 
the data analysis of this paper indicate that very 
few radio sources have a spectral slope index less 
than -1.5 between 20 and 92 cm. To ensure an 
NVSS detection for all sources with < —1.5 in 
a flux-limited sample requires iwENSS < 12-9 (see 



Eq. [7]). We positionally matched all £\venss < 
12.9 sources in the WENSS survey to NVSS and 
FIRST, and found that only 1% did not have a 
match in NVSS. We interpret this value as an up- 
per limit on the fraction of sources with a|o < 
— 1.5. Ninety-five per cent of the WENSS sources 
have both a FIRST detection and an NVSS de- 
tection within 30". The sample of NVSS-FIRST- 
WENSS matches with iwENSS < 12.9 in the uni- 
fied catalog contains 47,567 sources. Over 12,000 
of those have optical counterparts, from which a 
cosmological distance can be determined via pho- 
tometric redshift techniques. Photometric red- 
shifts from the SDSS are accurate to ~ 0.06 
for luminous red galaxies with redshifts z < 0.7 
(jPadmanabhan et al.l2005h . A further 3,000 of the 
radio-optical sample additionally have SDSS spec- 
troscopic redshifts. This dataset therefore repre- 
sents a major size increase compared to the sample 
of 327 sources used by Barai & Wiita. 



6. SUMMARY 

We have presented a unified catalog of radio ob- 
jects by combining observations from several ra- 
dio and optical surveys. The complete catalog 
is available online, along with scientifically use- 
ful data subsets. Using several parameters deter- 
mined from the measured radio fluxes, we investi- 
gated the distribution of sources in radio color- 
magnitude-morphology space and optical color- 
redshift-luminosity space. 

To create the catalog, we merged together 
the two 20 cm radio surveys, FIRST and NVSS. 
Included are FIRST-NVSS associations (three 
sources from each survey with the closest prox- 
imity to sources in the other survey, within 30"), 
FIRST sources with no NVSS counterpart, and 
NVSS sources without a nearby FIRST counter- 
part. Where sources had a counterpart in one 
of the other three surveys, we supplemented the 
catalog data with 6 cm, 92 cm, and optical obser- 
vations. 

For a preliminary data analysis, we converted 
the radio fluxes (6 cm, peak 20 cm, high- and 
low-resolution integrated 20 cm, and 92 cm) into 
several useful parameters describing 20 cm radio 
magnitude, 20 cm radio morphology, radio spec- 
tral slope, and optical identification/luminosity. 
The automatic morphology classification divides 
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sources into "complex", "resolved", and "com- 
pact" categories, with results that compare well to 
visual morphology classification. Using various ra- 
dio color-magnitude diagrams, we found that most 
radio sources have a steep radio spectral slope be- 
tween 6 and 20cm (a ~ —0.8, where F v oc v a ). 
However, the radio compact sources contain a sig- 
nificant fraction of flat-spectrum sources (a ~ 0). 
SDSS morphological identifications for optically- 
matched sources indicate that the steep-spectrum 
compact sources are primarily galaxies, while flat- 
spectrum compact sources are primarily optically 
unresolved (i.e., quasars). The complex and radio 
resolved sources contain comparable numbers of 
both galaxies and optical point sources. 

Focusing on sources with SDSS spectra, we in- 
vestigated the dependence of radio spectral slope 
and radio loudness on optical quantities (color, 
magnitude, redshift) for quasars and galaxies. For 
complex and radio resolved quasars, no signif- 
icant dependence of spectral slope on redshift 
or optical magnitude was found. For compact 
quasars, however, radio spectral slope seems to 
flatten out with increasing redshift/luminosity. 
The trend is consistent with a radio K-correction 
resulting from a spectrum that is steep at wave- 
lengths longer than 20 cm and flat for wavelengths 
shorter than 20 cm, but may also be an effect of 
orientation and Doppler beaming. Formally, a 
Kolmogorov-Smirnov test showed a dependence 
on optical color for the radio loudness distribution 
of steep-spectrum quasars. No color-dependence 
was found for the flat-spectrum quasars. In the 
galaxy sample, the compact morphology class 
tends to be bluer and fainter, with flatter radio 
slopes, than either the complex or radio resolved 
class. Among steep-radio-spectrum galaxies, the 
red sources have radio fluxes on average two orders 
of magnitude brighter than the blue sample. 

We have presented methods for selecting candi- 
date radio stars and candidate high-redshift galax- 
ies from the unified catalog. We concluded that a 
conservative color-selection attempting to reduce 
contamination from quasars still yields a sample 
that is completely dominated by quasars rather 
than radio stars. Spectra are therefore necessary 
to successfully identify radio stars. 

In a companion paper (A. Kimball et al. in 
preparation) , we will compare the results of several 
state-of-the-art models for radio galaxy evolution 



to observations yielded by the unified catalog, and 
test the unification paradigm for powerful radio 
galaxies and quasars. 

While this is not the first catalog to combine 
observations in the radio and the optical to fur- 
ther the study of radio galaxies and radio quasars, 
the unified catalog presented in this paper has sig- 
nificant advantages over previous studies. It is the 
largest catalog to date that combines optical ob- 
servations with multi-wavelength radio data. It 
can be used to define a long-wavelength (92 cm) 
flux-limited sample of nearly 50,000 sources de- 
tected in WENSS, FIRST, and NVSS. Comparing 
20 cm FIRST and NVSS fluxes allows us to place 
every one of those sources into one of the three ra- 
dio morphology classes, and the additional 92 cm 
detection allows us to calculate a radio spectral 
slope. Over one quarter of the long-wavelength 
flux-limited sample has an optical detection in the 
SDSS, including over 3,000 with spectra. The use- 
fulness of the catalog will increase as larger sur- 
veys at different wavelengths with deeper obser- 
vations become available. For example, the deep 
observations planned for the future Large Synop- 
tic Survey Telescope^ will provide a nearly 100% 
optical identification rate, and extremely accurate 
photometry in six bands. 

This material is based upon work supported un- 
der a National Science Foundation Graduate Re- 
search Fellowship, and by NSF grant AST-0507259 
to the University of Washington. The authors 
would also like to thank Wim de Vries and Paul 
Wiita for their helpful comments. 
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A. CATALOG MATCHING 



A.l. Matching FIRST and NVSS 

In this section, we describe in detail the process of matching sources in FIRST and NVSS, and present a 
method for selecting a clean sample of 20 cm sources. 

Before proceeding, it is necessary to define certain clarifying terms, some of which are data parameters 
found in the catalog. Our matching algorithm used a primary survey and a secondary survey. Source 
positions in the primary survey were used as search centers, and neighbors within a certain radius were 
found in the secondary survey. To find catalog sources we applied the matching algorithm twice: with 
FIRST as primary and NVSS as secondary, and again with NVSS as primary and FIRST as secondary. Each 
time, the three closest neighbors from the secondary survey were ranked by proximity to the primary source 
(within the pre-determined search radius). The advantage of performing the matching twice is that we can 
recover the three closest FIRST detections to a NVSS primary, and vice-versa. Each catalog row contains 
records for a single NVSS source and a single FIRST source. For an entry with an NVSS primary, the 
parameter matchflagjnvss is set to — 1, while matchflag_first is set either to 0, indicating a primary source 
with no FIRST match, or to 1, 2, or 3 indicating the proximity ranking of the FIRST match. Conversely, 
for FIRST primaries matchflag_first is set to— 1, and matchflagjnvss equals 0, 1, 2, or 3. Note that a single 
source can appear as the primary object in as many as three separate catalog entries. 

Due to the lower spatial resolution and higher flux limit of NVSS, individual NVSS sources can be 
resolved into separate components by FIRST. Therefore, multiple FIRST detections matched to the same 
NVSS primary are common, while multiple NVSS detections matched to the same FIRST primary are rare. 
Figure |2"21 gives a visual example of common types of objects found in the catalog. FIRST sources are shown 
as small squares and NVSS sources as large circles. In quadrant A, an unresolved NVSS source is resolved into 
three components in FIRST (example: a radio quasar with two lobes); a fourth FIRST detection lies nearby. 
When taking FIRST as the primary survey, each of the four FIRST sources would identify the NVSS object as 
the nearest neighbor. When taking NVSS as the primary survey, the three closest FIRST sources would match 
to the NVSS source while the fourth would not. The NVSS source in Quadrant A would appear in seven 
catalog entries. It would appear in four rows with matchflag_first = —1 and matchflagjnvss — 1. It would 
appear in three rows with matchflagjnvss = —1 and matchflag_first = 1, 2, or 3. To cull out the duplicates 
in order to select a clean sample, one would simply require matchflagjnvss = — 1 and matchflag_first = 1. 

The remaining objects in Figure [22] are simpler to match. Quadrant B contains a simple, isolated source 
identified by both surveys (examples: a radio quasar with no jets or a distant galaxy). Such an object would 
appear in two catalog rows: once with matchflag_first = —1 and matchflagjnvss = 1, and a second time 
with matchflagjnvss = —1 and matchflag_first = 1. Quadrant C contains a FIRST source without a NVSS 
counterpart (e.g., an object too faint to be detected by the latter); it would appear once in the catalog, with 
matchflag_first — —1 and matchflagjnvss = 0. In Quadrant D is a NVSS source unmatched by FIRST (e.g., 
a low surface-brightness galaxy). This source's single catalog entry would have matchflagjnvss — —1 and 
matchflag_first = 0. 

To find all sources detected by both surveys while minimizing duplicates, the simplest selection is to 
require matchflagjnvss = —1 and matchflag_first = \. The resulting data set consists of 142,622 sources in 
the overlap region, with only 158 (0.11%) duplicate NVSS sources matched to the same FIRST source. 
Specific examples of data selections are given in Appendix [B] 

A.2. Matching to the SDSS 

To match to the SDSS photometric catalog, we implemented a different technique from that used in 
matching the radio surveys. The motivation was to find the nearest neighbor, as well as the brightest 
neighbor within a pre-defined distance D, where D was adjusted according to the distance of the nearest 
neighbor. We began by finding the nearest optical neighbor within 60" of each catalog entry. Where the 
distance d near to the nearest neighbor was less than 3", we defined D — 3". For 3" < d near < 10", we 
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defined D = 10", for 10" < d near < 30" we denned D = 30", and for 30" < d near < 60" we denned D = 60". 
We then recorded the brightest neighbor within distance D. Each catalog entry contains two sets of data 
parameters corresponding to SDSS photometric matches: one set is identified by the prefix "near" and the 
other by the prefix "bright" (e.g. ll near_modelmag_u" , "bright jmodelmag-u") . Thirty-five per cent of FIRST 
sources in the overlap region have an SDSS counterpart within 3", 61% within 10", and over 98% within 
30". In 98.8% of cases where an optical match was found within 3", the brightest and the nearest SDSS 
objects within 3" are the same. 

B. SELECTING DATA SAMPLES 

In this section, we define in detail the samples described in this paper, and present several scientifically- 
useful data subsets available for download on the catalog website0 Samples described in this appendix but 
not available online can be assembled from the complete catalog using the selection requirements given here. 

Table [8] presents a summary of the data sets and their selection requirements. The meanings of catalog 
parameters referenced in the table (in italics) are as follows. Parameters ra and dec are right ascension 
and declination, respectively, in decimal degrees. Overlap indicates a source inside (value 1) or outside 
(value 0) the catalog overlap region. Matchflagjnvss and matchflag_first are set to —1 to indicate primary 
source catalog, to indicate an isolated source, or 1, 2, 3 to rank matched source by proximity to primary 
source (see Appendix [A~1 for details). Distance indicates separation between FIRST and NVSS neighbors, in 
arcseconds. First_id, wenssJd and gb6J,d are unique identifiers for the FIRST, WENSS and GB6 sources, 
respectively. W ens s -distance and gb6-distance indicate separation, in arcseconds, between WENSS or GB6 
neighbors and the 20 cm source. Near_distance and nearJ,ype indicate respectively the distance to and object 
type of the nearest neighbor in the SDSS photometric survey, while spec-type indicates the SDSS spectral 
classification for objects from the spectroscopic survey. MatchtotSO indicates the number of neighbors within 
30" of the primary 20 cm survey source (where neighbors come from the secondary 20 cm survey). For a full 
description of all 117 catalog parameters, we refer the reader to the catalog website. 

B.l. Data sets discussed in the paper 

The data samples discussed in this paper consist of sources detected by both FIRST and NVSS, sources 
detected by the three largest radio surveys (FIRST, NVSS, WENSS), sources detected by all four radio 
surveys (FIRST, NVSS, WENSS, GB6), and sources detected by all four radio surveys as well as the SDSS. 
Spectroscopic galaxy and quasar samples are created from the FIRST-NVSS- WENSS sources. 

Sources detected by both FIRST and NVSS, matched within 25", are listed as sample C in Table [8] An 
efficient way of finding these pairs while culling duplicates is to select NVSS sources, and their nearest FIRST 
neighbor. These sources can be selected by requiring overlap — 1, matchflagjnvss =— 1, matchflag_first = 1, 
and distance < 25. The resulting sample contains 141,881 objects (48deg -2 ) detected by both 20cm surveys, 
which can therefore be classified by radio morphology. 

Sample D contains sources detected by FIRST, NVSS, and WENSS. This is the subset of sample C with 
a WENSS counterpart within 30". Additional requirements to select this data set are wenssjid^ and 
wenss_distance < 30. The resulting sample contains 63,660 sources (21.5 deg~ 2 ). These sources have three 
20 cm fluxes and a 92 cm flux measurement, and can therefore be categorized by both radio morphology and 
spectral slope (afo)- 

Sample E contains sources detected by all four radio surveys; it is the subset of sample D that has a 
GB6 counterpart within 70". Selecting these sources requires additional parameter settings gb6-id^ and 
gb6-distance < 70. The sample is comprised of 12,414 sources (4.2 deg -2 ), which can be categorized by radio 
morphology and two spectral slope measures. 



11 http://www.astro.washington.edu/akimball/radiocat/ 
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Sample F contains sources detected by all four radio surveys and the SDSS. To select these sources from 
sample E requires near_type ^ and near_distance < 2. The sample contains 4,732 sources (1.6 deg -2 ), which 
can be categorized by radio morphology and two spectral slope measures, and identified as galaxies or point 
sources. 

Sample G contains spectroscopically-identified galaxies detected by NVSS, FIRST, WENSS, and SDSS. 
These sources are a subset of Sample D, requiring spec-type = 2 and near_distance < 2. It contains 2,885 
galaxies. Note that in §4.4.2[ we restricted our analysis to a smaller, flux- limited sample of 1,656 spectroscopic 
galaxies with r < 17.77. 

Defining the spectroscopic quasar sample was more complicated than defin ing the spectroscopic galaxy 
sample. As discussed in ^2.2.21 the quasar catalog of ISchneider et alj (|2007l) is cleaner and hence more 
reliable than the set of spectroscopically-identified quasars supp lied by SDSS DR6. To identify radio quasars 
with spectra, we therefore opted to match radio sources to the ISchneider et al.l <|2007l ) catalog, temporarily 
ignoring the spectroscopic identifications obtained from SDSS DR6. The sample we analyzed (sample H) 
consists of 1,288 NVSS-FIRST- WENSS sources with a spectroscopic quasar within 2". This file is available 
only in ASCII format with 105 columns, and contains data not available in the full catalog. The first 74 
columns contain data specific to the DR5 known quasar catalog; the remaining 31 rows contain radio fluxes 
and other data parameters from the unified catalog. A header in the file gives a description of all 105 
columns. 



B.2. Other data sets available online 

Other data files available online provide the complete catalog, a small spatial subset containing all data 
parameters, isolated FIRST-NVSS sources, isolated FIRST-NVSS sources with a counterpart in WENSS or 
GB6, and isolated FIRST-NVSS sources with optical counterparts. 

The complete catalog (sample A) is available as a compressed tar archive. It contains 72 files (in fits 
format), each covering a strip of sky with a width of 5°in right ascension. This archive contains all 117 data 
parameters for each of the 2,724,343 catalog entries; no selections of any kind were applied to this sample. 
A full explanation of all the data parameters can be found on the website. 

A small subset (sample B) of the catalog is also available, intended for testing purposes, to demonstrate 
the catalog format. It covers approximately 106 deg 2 of sky in the range 150° < R. A. < 165° and 40° < dec. < 
50°. The file contains all 117 data parameters for the 16,453 catalog entries in that sky region. 

We provide a set of isolated FIRST-NVSS matches within the catalog overlap region (sample I). The set 
of isolated sources is comprised of 109,825 NVSS-FIRST detections within a matching radius of 15", and no 
other FIRST neighbors within 30". The file includes positions (right ascension and declination of the FIRST 
source), distance between the FIRST and NVSS positions, and 20cm fluxes. For sources with a match in 
GB6 or WENSS within 120" (47,490 total), it also contains distance to the matched source along with its 
flux. Sources without a match have the corresponding data parameters equal to zero. 

We provide the subset of isolated FIRST-NVSS matches with a SDSS counterpart within 2" (sample J). 
There are 44,851 such sources in the catalog overlap region. In addition to data parameters describing FIRST 
and NVSS distance and fluxes, this file contains the position, photometric type, magnitudes, and magnitude 
errors for the nearest SDSS match. 

Finally, we provide the set of candidate 9,619 high-redshift galaxies discussed in §5.21 (sample K). 
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Fig. 1. — Using sparse sampling, we show the sky 
coverage of each of the surveys used to created 
the matched catalog. The blue line shows the 
2955 deg 2 region of overlap of the radio surveys 
and the SDSS photometric survey. The 2894 deg 2 
region outlined by the orange and blue lines shows 
the region included in the spectroscopic sample 
of SDSS DR6. The dashed red line indicates the 
Galactic plane. 
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Fig. 2. — Differences in equatorial coordinate po- 
sitions (top: right ascension; bottom: declina- 
tion) for FIRST-SDSS and NVSS-SDSS matches. 
The full-width at half-maximum (FWHM) of the 
FIRST-SDSS distribution is about 1". For the 
NVSS-SDSS distribution the FWHM is about 6". 
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Fig. 3. — Diamonds indicate distributions of 
distance between FIRST and (A) NVSS, (B) 
WENSS, (C) GB6, and (D) SDSS positions 
for close pairs from sub-regions of the catalog. 
Squares show the random match distribution, es- 
timated by off-setting the NVSS, WENSS, GB6, 
or SDSS positions by 1° in right ascension. Panel 
insets show the completeness (solid curve) and ef- 
ficiency (dotted curve) as a function of matching 
radius, estimated by a model fit to the nearest 
neighbor distribution. The solid vertical lines indi- 
cate radii used to define catalog matches. Dashed 
vertical lines show the more conservative matching 
radii used for the analysis presented in this paper. 
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Fig. 4. — Sky density of FIRST sources matched 
to the other radio surveys as a function of magni- 
tude. The total sky density of each sample is indi- 
cated below the legend in each panel. This plot in- 
cludes FIRST components from the same complex 
source; i.e. multiple FIRST detections matched to 
a single NVSS source. Panel A: FIRST sources 
that were or were not detected in NVSS. Sky den- 
sities for sources brighter than lOmJy (13.9 mag) 
are 18deg~ 2 (triangles) and 0.53 deg~ 2 (squares). 
Panel B: FIRST-NVSS sources with and without 
a counterpart in GB6. Sky densities for sources 
brighter than 100 mjy (11.4 mag) are 1.5 deg -2 
(triangles) and 0.43 deg -2 (squares). Panel C: 
FIRST-NVSS sources with and without a counter- 
part in WENSS. Sky densities for sources brighter 
than lOmJy (13.9 mag) are 16 deg -2 (triangles) 
and 1.7deg -2 (squares). 
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Fig. 5. — Top: Correlation between NVSS (lack- 
ing FIRST counterpart within 30") and WENSS 
sources (diamonds); random match distribution 
(squares) is shown for comparison, measured by 
offsetting WENSS positions by 1° in right as- 
cension. The inset shows completeness (solid 
curve) and efficiency (dotted curve) as a func- 
tion of matching radius. Solid vertical lines in- 
dicate matching radius used to create the cata- 
log; dashed vertical lines show the matching radius 
used for this paper's analysis. Bottom: Sky den- 
sity of NVSS sources without FIRST counterpart 
as a function of magnitude. The total sky den- 
sity is indicated below the legend. Sky densities 
for sources brighter than lOmJy (13.9 mag) are 
0.35 deg~ 2 (triangles) and 0.21 deg -2 (squares). 
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Fig. 6. — Top: Difference in 20 cm radio magni- 
tudes measured by FIRST and NVSS as a function 
of iNVSS- Where individual points saturate the 
plot, we show the distribution using contours. The 
dashed horizontal lines indicate peaks at At = 
and At = 0.7. The solid horizontal line shows 
the At — 0.35 separator between between complex 
and simple sources. Vertical dashed lines mark the 
magnitude bins shown in the lower panel. Bot- 
tom: Distribution in At for four tNVSS bins. The 
vertical solid line shows the At = 0.35 separator 
between complex and simple sources. 
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Fig. 7. — Top: Distribution of At for sources with 
fcrvss < 13.5. The dashed line indicates the sepa- 
ration between the simple and complex morphol- 
ogy classes. Middle: Distribution of 0, an esti- 
mate of source size. The solid curve corresponds 
to simple sources; the dotted curve corresponds 
to complex sources. The vertical line separates 
unresolved and resolved sources at 20 cm. Bot- 
tom: Distribution of At vs. Iog0 2 for sources 
with £nvss < 13.5. The dashed lines indicate 
the adopted separation of morphology classes: the 
complex sample is above the horizontal dashed 
line, the compact sample is left of the vertical 
dashed line, and the resolved sample is to the right 
of the vertical dashed line. 
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Fig. 8. — Examples of radio sources classified by 
morphology (FIRST 1' x 1' stamps normalized 
to the brightest pixel, shown with square root 
stretch). The top row contains compact radio 
sources, the second row shows resolved sources, 
and the remainder are complex. The left two 
columns show sources optically identified as galax- 
ies by the SDSS; the right two columns show 
SDSS-identified quasars. 




Fig. 9.— Top: Distribution of NVSS-FIRST 
sources classified according to radio morphology, 
as a function of £nvss- Bottom: Fraction of 
sources brighter than t max belonging to each mor- 
phology class. 
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Fig. 10. — Color-magnitude (left), color-color 
(middle), and median magnitude (right) diagrams 
for sources detected in all four radio surveys and 
classified by radio morphology (top: complex; 
middle: resolved; bottom: compact). Vertical 
dashed lines show the o^ 2 = —0.5 separation be- 
tween steep- and flat-spectrum sources in the 20- 
92 cm range. Left column: The distribution of 

20- Center 
The 

diagonal, dashed red lines show the a^o = a 6° 
locus. Right column: Color-color diagrams with 
bins color-coded according to the median value 
of NVSS magnitude t^vss- The strong gradi- 
ents in the right column are due to steep-spectrum 
sources, which are biased to bright sources by the 
shallow GB6 survey. 



NVSS magnitude vs. spectral slope a^Q 
column: af vs. q^q color-color diagrams 
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Fig. 11. — The magnitude distribution of compact 
morphology radio sources, observed by all four ra- 
dio surveys and brighter than £nvss = 12, catego- 
rized into four subclasses by spectral shape. Top: 
Representations of a standard spectral shape in 
each class, defined by radio fluxes at 6, 20, and 
92 cm. Steep and flat refer to sources with mono- 
tonic radio flux measurements. Steep sources have 
ag 2 < —0.5; flat sources have ag 2 > —0.5. Peaked 
refers to sources which are brightest at 20 cm. In- 
verted refers to sources which are faintest at 20 cm. 
Middle: Twenty centimeter magnitude distribu- 
tion for each spectral shape subclass. Bottom: 
Fraction of sources brighter than f max belonging 
to each spectral shape class. 
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Fig. 12. — Magnitude distributions for sources de- 
tected by all four radio surveys. Panel A: Full 
sample divided by radio morphology. Panel B: 
Full sample divided by optical identification. 
Panel C: Optically-identified subset divided by ra- 
dio morphology. Panels D-F: Complex, resolved, 
and compact radio sources, respectively. 
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Fig. 13. — Radio color-color diagrams for sources 
detected by all four radio surveys divided into the 
nine radio/optical classes. The abcissa is a^oi the 
ordinate is a|°. From top to bottom, rows con- 
tain complex, resolved, and compact radio sources. 
From left to right the columns contain optically 
undetected sources, optically resolved sources, and 
optically unresolved sources. Contours in each 
column show the underlying distribution of all 
sources in the corresponding column. Contour lev- 
els and histogram normalization are arbitrary for 
each column. 

#/deg7bin (arbitrary normalization for each column) 
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Fig. 14. — Distribution of the average of spec- 
tral indices a§° and a^o for £nvss < 12 sources 
detected in all four radio surveys. The vertical 
dashed lines show the separation between steep- 
and flat-spectrum sources. Note that the Y-axis 
range varies. Panel A: Full sample divided by ra- 
dio morphology. Panel B: Full sample divided 
by optical identification. Panel C: Optically- 
identified sources divided by radio morphology. 
Panels D-F: Complex, resolved, and compact ra- 
dio sources, respectively. 
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Fig. 15. — Radio color-magnitude diagrams for 
the nine radio/optical subclasses from the NVSS- 
FIRST-WENSS sample, categorized by radio mor- 
phology and optical identification. From top to 
bottom, the rows contain complex, resolved, and 
compact radio sources. From left to right, the 
columns contain the 46,517 optically undetected, 
the 12,780 optically resolved, and the 4,363 op- 
tically unresolved sources. The left axis is la- 
beled with radio magnitude; the right axis with 
corresponding flux. Contours correspond to all 
the sources in each column; contour levels and 
histogram normalization are arbitrary for each 
column. Note that radio sources without an 
SDSS match are dominated by complex steep 
sources (a^o < —0.5), while those matched to an 
SDSS point source are dominated by compact flat 
sources (o^o > —0.5). 
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Fig. 16. — Magnitude distributions for sources de- 
tected by NVSS, FIRST, and WENSS. Panel A: 
Full sample divided by radio morphology. 
Panel B: Full sample divided by optical identifica- 
tion. Panel C: Optically-identified sources divided 
by radio morphology. Panels D-F: Complex, re- 
solved, and compact radio sources, respectively, 
divided by optical identification. 
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Fig. 17. — Distributions in absolute magnitude- 
redshift-spectral index parameter space for the 
spectroscopic quasar sample. The dotted and 
dashed lines indicate the i < 20.2 and i < 19.1 
limits, respectively, used in the quasar target selec- 
tion algorithm. Panel A: absolute magnitude vs. 
redshift for the full sample, with bins in redshift 
and apparent magnitude color-coded by the me- 
dian value of ago (upper colorbar). The histogram 
at the bottom of the panel shows the redshift dis- 
tribution. Panels B through D: panel A divided 
into the three radio morphology subclasses. The 
percentage of quasars in each morphology class is 
given. Full histograms show the redshift distri- 
bution in each panel; for comparision the redshift 
distribution of the full quasar sample is also shown 
(dotted line). Histograms have been normalized 
such that all have equal area. Panels E and F: 
Spectral index vs. optical absolute magnitude dis- 
tributions for complex and compact radio sources, 
respectively. Points are color-coded by redshift. 
Circles show the median value of a^o in bins of 
1 magnitude. 
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Fig. 18. — A simulated radio spectrum for which 
a|o would appear flat at high redshift (z) and 
steep at low redshift. Observed 20 and 92 cm 
fluxes correspond to shorter rest-frame wave- 
lengths for a high-redshift object. 
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Fig. 19. — Distribution of af 2 , a measure of ra- 
dio loudness, for the spectroscopic quasar sam- 
ple. The top two panels show the region of ra- 
dio magnitude-optical magnitude space selected 
so as not to bias toward radio quiet or radio 
loud sources. The lower two panels show a? 2 his- 
tograms for red (u — r > 0.7; solid line) and blue 
(u — r < 0.7; dashed line) quasars. The left pan- 
els correspond to sources with steep radio spectra; 
panels on the right correspond to sources with flat 
radio spectra. 
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Fig. 20. — Distribution in color, magnitude, spec- 
tral index parameter space for the flux-limited 
(r < 17.77) spectroscopic galaxy sample. Panel A: 
Color-magnitude diagram for the full sample with 
bins color-coded by the median value of a^p . Con- 
tours show the distribution of all r < 17.77 
SDSS spectroscopic galaxies. The histogram be- 
low shows the u — r distribution. Panels B 
through D: panel A divided into the three radio 
morphology classes, as labeled. Points are color- 
coded by Qf^o- Percentage of galaxies in the sub- 
class is given in the lower right corner. Dashed 
lines at u — r = 3 and M r = —22.5 are provided to 
guide the eye. The histogram underneath shows 
the u — r distribution (solid line), compared to 
the distribution for the full sample (dotted line; 
equivalent to solid line histogram in Panel A). 
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Fig. 21. — Distributions in color, magnitude, and 
radio to optical spectral slope for the flux-limited 
(r < 17.77) spectroscopic galaxy sample. Pan- 
els A and B: Color-magnitude diagrams (M r vs. 
u — r) for steep- and flat-spectrum radio sources, 
respectively; points are color-coded by the radio 
to optical spectral slope (a? 2 ). The solid line his- 
togram shows the u — r distribution; the dotted 
line histogram shows the u — r distribution of the 
full sample (equivalent to solid line histogram in 
Panel A of Figure [20]) . The dashed vertical line 
indicates u — r = 2.5. The percentage of galax- 
ies in each class is given in the lower right corner. 
Panels C and D: Normalized distributions of a\ 2 
for red (u — r > 2.5) and blue (u — r < 2.5) radio 
galaxies from the top two panels. 
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Fig. 22. — Examples of different types of sources 
found by matching NVSS and FIRST. Circles rep- 
resent NVSS detections; squares represent FIRST 
detections. The sizes of the two types of symbols 
help illustrate the different spatial resolution of 
the two surveys. Quadrant A contains a complex 
source resolved into three components by FIRST 
but unresolved by NVSS, as well as a fourth (un- 
related) FIRST detection. Quadrant B contains a 
simple source identified by both surveys, with no 
nearby neighbors. Quadrant C contains a FIRST 
detection with no NVSS counterpart, while Quad- 
rant D contains a NVSS detection with no FIRST 
counterpart. 
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Table 1 

Technical comparison of the radio source surveys. 



statistics 


FIRST 


NVSS 


GB6 


WENSS 


Frequency [MHz] 


1400 


1400 


5000 


326 


Wavelength [cm] 


20 


20 


6 


92 


Resolution [arcsec] 


5 


45 


210 


54 


Sky coverage 


b > 30° 


5 > -40° 


< S < 75 


5 > 30° 


Sky area [deg~ ] 


10,000 


33,000 


17,000 


10,000 


Lim. Flux density 


1 


2.5 


18 


18 


(5<7 rms , mJy) 










Lim. AB mag 


16.4 


15.4 


13.3 


13.3 


Source density [deg~ 2 ] 


97 


55 


4.4 


22 



Table 2 

Completeness and efficiency as a function of matching radius 



Matched surveys 


Radius [arcsec] 


Completeness 


Efficiency 


FIRST-NVSS 


25 


0.990 


0.96 




30 


0.997 


0.96 


FIRST-WENSS 


30 


0.990 


0.92 




120 


1.000 


0.74 


FIRST-GB6 


70 


0.979 


0.79 




120 


1.000 


0.68 


FIRST-SDSS 


2 


0.982 


0.95 




3 


0.998 


0.90 




10 


1.000 


0.47 


NVSS- WENSS 


55 


0.910 


0.85 




120 


0.995 


0.57 



Note. — To calculate these values, the nearest neighbor distribu- 
tions were fit to a Gaussian (representing physical matches) plus a 
rising linear function (representing random matches). 
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Table 3 

Double-Gaussian fits to FIRST-NVSS At distributions 



Bin Simple sources Complex sources 





Fraction a 


Mean b 


CT C 


Fraction a 


Mean b 


CT C 


^NVSS < 12 


0.59 


0.035 


0.044 


0.41 


0.58 


0.35 


12 < i NVS s < 13 


0.67 


0.064 


0.074 


0.33 


0.64 


0.37 


13 < i N vss < 14 


0.65 


0.12 


0.13 


0.35 


0.75 


0.47 


14 < i N vss < 15 


0.59 


0.23 


0.24 


0.41 


0.75 


0.49 



Note. — The distributions are shown in the lower panel of Figure [5] 
a Fraction of total area belonging to each Gaussian. 
b Position of Gaussian mean in At. 
c Standard deviation. 



Table 4 

Comparison of automatic and visual morphology classification 



Visual 


Total 


Complex (%) a 


Resolved (%) a 


Compact (%) a 


Classification 


QSOs 


Gals. 


QSOs 


Gals. 


QSOs Gals. 


QSOs 


Gals. 


Complex 


202 


286 


77 


95 


8 4 


15 


1 


Resolved 


87 


79 


8 


42 


72 58 


20 





Compact 


273 


144 


5 


26 


9 26 


80 


48 



Note. — Results compiled for a set of 1071 radio quasars and radio galaxies (identified 
spectroscopically from the SDSS). Visual classification is listed at the left; automatic 
classification is listed along the top. 

a Numbers refer to the percentage of quasars/galaxies in each visual classification cate- 
gory which received a particular automatic classification. 
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Table 5 

Gaussian fits to (a^l + al°)/2 distributions 



Subclass deg 2 mean a median FWHM a 



Complex 13 


















All 


1 


.3 


-0. 


.868 


-0 


.863 





.433 


no SDSS 


1 


.1 


-0. 


.886 


-Q 


.881 





.406 


galaxy 


0. 


.15 


-0. 


.742 


-0 


.732 





.433 


pt. source d 


0. 


.069 


-0. 


.705 


-0 


.662 





.629 


Resolved 13 


















All 





.08 


-0. 


.850 


-0 


.843 





.433 


no SDSS 





.38 


-0. 


.895 


-Q 


.892 





.416 


galaxy c 


0. 


.19 


-0. 


.808 


-0, 


.805 





.385 


pt. source d 


0. 


.11 


-0. 


.742 


-0 


.733 





.571 


Compact 13 


















All 





.99 


-0. 


,713 


-0 


.680 





.822 


no SDSS 


0. 


.45 


-0. 


,839 


-0 


.815 





.530 


galaxy c 


0. 


.24 


-0. 


.710 


-0, 


,670 





.682 


pt. source d 


0. 


.30 


-0. 


,302 


-0, 


.288 


1.03 



Note. — Values correspond to histograms shown 
in Figure [T4l Sources were detected by all four radio 
surveys, and follow the additional constraint £nvss < 
12. 

a Estimated from a Gaussian fit to the (o^q+chq ) /2 
distribution. 

b 20cm radio morphology. 

c Unresolved by SDSS. 

d Resolved by SDSS. 
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Table 6 

Gaussian fits to a™ distributions 



Subclass 


deg" 2 


mean a 


median 


FWHM a 


Complex 13 














All 


1.3 


-0. 


.786 


-0 


.782 


0.323 


no SDSS 


1.1 


-0. 


.798 


-Q 


.797 


0.301 


galaxy 


0.15 


-0. 


.659 


-0 


.654 


0.323 


pt. source d 


0.069 


-0. 


.718 


-0 


.690 


0.396 


Resolved" 3 














All 


0.68 


-0. 


,764 


-0 


.765 


0.339 


no SDSS 


0.38 


-0. 


,807 


-0 


.817 


0.320 


galaxy c 


0.19 


-0. 


,700 


-0, 


.702 


0.311 


pt. source d 


0.11 


-0. 


,713 


-0, 


,710 


0.378 


Compact 13 














All 


0.99 


-0. 


,632 


-0 


.563 


0.640 


no SDSS 


0.45 


-0. 


,727 


-0 


.683 


0.440 


galaxy 


0.24 


-0. 


,593 


-0, 


.553 


0.576 


pt. source d 


0.30 


-0. 


,271 


-0, 


.229 


1.05 



Note. — Values correspond to the data samples in- 
vestigated in Fig. Q3] and Table Sources were de- 
tected by all four radio surveys, and follow the ad- 
ditional constraint £nvss < 12- The artjci histograms 
are not shown. 

a Estimated from a Gaussian fit to the distribution. 

b 20 cm radio morphology. 

c Unresolved by SDSS. 

d Resolved by SDSS. 
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Table 7 

Candidate radio star selection 



Data Source 


All a 


in FIRST b 


i < 19.0 C 


Non-QSO colors' 1 


w/ spec. 6 


stellar spec. f 


opt. point sources 6 


19,133,672 


12,100 


2,424 


532 


406 


78 


quasar catalog 6 


35,450 


2,927 


1,761 


289 


289 


2 




(37,010) 


(3,147) 


(1,931) 


(358) 


(358) 


(2) 



Note.— 

a All sources in the catalog overlap region. (For the optical point source sample, there are 5,841,685 
sources with i < 19). 

b Within 1 arcsec from a FIRST source. 

c Using model magnitudes (1,653 for quasar catalog when using psf magnitudes.) 
d Selected by rejecting sources with u — g < 0.8 AND —0.2 < g — r < 0.6. 
6 Selected from SDSS Data Release 6. 

f Visually-inspected to confirm automatic classification as stellar spectrum. 

g Schncidcr et al. (2006) catalog corresponding to an unresolved SDSS source (values in parentheses 
correspond to the entire quasar catalog). 
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Table 8 

Catalog subsets available for download 



Sample Description Selection rows size 3. 



A Complete catalog 
B 100 deg 2 subset 



C Detected by FIRST 
and NVSS 



None 

140° <ra< 165° 
40° <dec< 50° 

overlap = 1 
matchflag-nvss = — 1 
matchflag-first = 1 
distance < 25" 



2,724,343 257M (1.4G) 
16,453 3.4M (15M) 

141,881 5.4M (8.7M) 



D Detected by FIRST, 
NVSS, and WENSS 



(subset of C) 
wenss-flux ^ 
wenss-distance < 30" 



63,660 



E Detected by FIRST, 
NVSS, WENSS, 
and GB6 

F Detected by FIRST, 
NVSS, WENSS, GB6, 
and SDSS 

G Spectroscopic galaxies 
detected by FIRST, 
NVSS, WENSS, 
and SDSS 



(subset of D) 
gb6-flux ^ 
gb6-distance < 70" 

(subset of E) b 
near-type ^ 
near-distance < 2" 

(subset of D) b 
spec-type = 2 
near-distance < 2" 



12,414 



4,732 



2,885 



625K 



439K 



H Spectroscopic quasars 
detected by FIRST, 
NVSS, WENSS, SDSS 

I Isolated d FIRST-NVSS 



Isolated d FIRST-NVSS, 
and SDSS 



(subset of C) c 
Match to Schneider et 
al. (2007) within 2" 

(subset of C) 
distance < 15" 
No other FIRST 
sources within 30" 

(subset of I) b 
near-type ^ 
near-distance < 2" 



1,288 



755K 



109,825 4.1M (6.8M) 



44,851 3.7M (5.7M) 
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Table 8 — Continued 



Sample Description 



Selection 



# rows 



K Candidate high- 
redshift galaxies 

(*>1) 



a 



(subset of E) 
(near_ra=0 OR 
near_distance > 3") 
11 < -0.5, af < -0.5 
remove suspicious 
objects (gSU 



9,619 2.5M (5.1M) 



Note. — The first column gives the subset identifier, the second column con- 
tains a brief description, the third column lists the requirements to select the 
data set from the complete catalog, and the fourth column shows the number of 
sources in the data set. The fifth column indicates the size of the downloadable 
file for datasets available (with limited data parameters) on the catalog website. 
For compressed files, the uncompressed size is listed in parentheses. 

a File size. For compressed files, number in parentheses refers to uncompressed 
size. For samples which are subsets of downloadable files (in terms of both selected 
objects and relevant data parameters), no size is given. 

b While these objects are a subset of another sample, the file includes additional 
SDSS data parameters. 

c The sample of spectroscopic quasars is a subset of ob jects i n Sam ple C; however 



the file includes data parameters from the ISchneider et alJ (|2007T ) catalog which 
are not part of the unified radio catalog. 

isolated in this context refers to FIRST-NVSS pairs with no other FIRST 
sources within 30". 
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